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^ Abstract 

Q We study the low-temperature properties of a spin-^ magnetic impurity cou- 

pled to a one-dimensional interacting electron system. Using the newly de- 
veloped formalism by Affleck and Ludwig, with a scale invariant boundary 
5_j ■ condition replacing the impurity, we exploit boundary conformal field the- 

ory to deduce the impurity thermal and magnetic response. In the case of 
only forward electron scattering off the impurity, we predict the same criti- 
cal scaling as for the two-channel Kondo effect for non-interacting electrons, 
but with a novel Wilson ratio. Backward electron scattering off the impurity 
destabilizes this behavior and drives the system to a new fixed point. In the 
case of equal amplitudes for forward- and backward scattering (Kondo inter- 
action), we show that there are only two types of scaling behaviors consistent 
with the symmetries of the problem: either a local Fermi liquid or a critical 
theory with an anomalous specific heat. The latter case agrees with a recent 
"poor-man-scaling" result proposed by Furusaki and Nagaosa. 
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I. INTRODUCTION 



Quantum many-particle systems sometimes exhibit a growth of an effective coupling at 
low energies, resulting in a non-perturbative ground state. The maybe simplest example of 
this phenomenon is the Kondo effect J5J, arising from the exchange interaction between a 
spin-| magnetic impurity and a gas of free quasi-particles ("dressed electrons") in an s-wave 
band. As the temperature is decreased, the system crosses over from weak electron-impurity 
coupling to strong coupling, with a complete screening of the impurity spin at T = 0. The 
resulting ground state is of Fermi liquid type, with the single quasi-particle wave functions 
acquiring a phase shift ("one-channel Kondo effect") 0. The picture changes when electrons 
in degenerate orbital bands are allowed to interact with the impurity, and the ground state 
is now described by a non- Fermi liquid fixed point ( "multi-channel Kondo effect" ) |3| . 

What is the corresponding scenario for an interacting one-dimensional electron system 
coupled to a magnetic impurity? The question may soon become of experimental relevance, 
considering the rapid progress in the fabrication and study of very narrow conduction chan- 
nels ("quantum wires"), obtained, for example, by gating 2D electron gases in GaAs inversion 
layers Q. A possible laboratory realization would be a single conduction channel with a 
trapped atom containing two (or several) spin levels. (The related problem of tunneling 
through potential barriers in ID correlated systems is already being addressed by experi- 
mentalists ||, following the pioneering work of Kane and Fisher 0.) The question is also 
interesting considering recent work on exotic superconductivity, where the analogue to the 
multi-channel Kondo effect has been exploited |7j. Other realizations of Kondo physics that 
have been proposed include two- level tunneling in metallic glasses |§, and certain heavy- 
fermion materials ||. Treating the effect of a magnetic impurity in the presence of interacting 
electrons may offer a new perspective on these intriguing connections. 

More importantly, one is here faced with an archetype problem of describing the interplay 
between direct fermion correlations (from interaction and statistics) and correlations induced 
via a coupling to a local quantum mechanical degree of freedom. As is well known, the 
notion of free quasi-particles (low-temperature Fermi liquid) breaks down in one dimension: 
any arbitrarily small electron-electron interaction wipes out the single-particle poles of the 
electron propagator, leaving behind only collective charge- and spin-density excitations. In 
the limit of weak electron-electron interaction, these excitations are well described by a 
system of non-interacting spin-charge separated bosonic modes (Luttinger liquid) |To| , pT| . 
By adding a localized magnetic impurity, one confronts the problem of how to incorporate 
its coupling to single electrons in the description of the bosonic collective degrees of freedom. 

A first attack on the problem was launched by Lee and Toner [ 12]] , employing Abelian 
bosonization followed by a perturbative scaling analysis of a resulting "kink-gas" action. 
For the case of a spin-| impurity, and with the electron gas away from half-filling, it was 
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found that the Kondo temperature Tk (setting the scale for the weak-to-strong coupling 
crossover) depends on the bare Kondo coupling A^- in a power-law fashion, Tk ~ {^kTo) 2 ■ 
Here r is a short-time cutoff, and rj is the exponent characterizing the equal-time spin- 
spin correlation in a Luttinger liquid. For temperatures T < T K , the physics is controlled 
by some strong-coupling fixed point — as in the ordinary Kondo problem — not directly 
accessible via this kind of analysis. In a recent work, Furusaki and Nagaosa [I3~| derived a 
set of improved scaling equations in the weak-coupling regime, preserving the spin SU(2) 
symmetry of the problem. By tentatively extending these equations to the strong-coupling 
regime, Furusaki and Nagaosa conjectured that the fixed-point Hamiltonian consists of two 
semi-infinite Luttinger liquids and a completely screened impurity (decoupled spin singlet). 
The low-temperature impurity contributions to the specific heat and magnetic susceptibility 
were calculated to C imp ~ T( 1 ' Kp '~ 1 and Ximp ~ T°, respectively, with K p the usual Luttinger 
liquid charge parameter [TOLpl. Support for this scenario can be found in earlier work p4[ on 



impurity spins in antiferromagnetic spin- chains (a "stripped-down" version of the Kondo 
effect in a Luttinger liquid). In the case of an external S — | impurity coupled to a single 
site on the chain, the impurity was found to be completely screened, severing the chain at 
the impurity site. 

In this paper we explore the problem using exact methods, expanding upon results an- 
nounced in |15[]. We shall begin by studying a simplified model. Specifically, we consider a 



spin- 1 impurity coupled with equal strength to two next-nearest neighbor sites on a Hub- 
bard chain. In the continuum limit, and with quarter-filling of the band, this becomes a 



Tomonaga-Luttinger (TL) model (|T6| , |T7i with forward electron-impurity scattering only: 



7i = Ti-TL + 7~L 



(1-1) 



where 



and 



H 



TL 



1 f* d d 

dx { v F [. ^ Lta {x)i—iP L (T {x) : - : ^M) 1 ^ r,A x ) 



2tt 



dx 



;s=L,R 



+ 9 ■^R,ai X )i ) L,A X )i ) L,-A X )i ) R,-a{ X )^ 



;i-2) 



7Yp = A 



1.3) 



Here ip La {x) and ip Ra (x) are the left- and right-moving components, respectively, of the 
electron field ^ a (x), with spin projection a =f , |, expanded about the Fermi momenta ±kp 
in the long- wavelength limit: 
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= e-^ L » + e lk ^ R Jx). (1.4) 
The fields are normalized such that 

{VV,<r( a 0>V'L(2/)} = 2ir6 r8 6 t , tl 6{x -y) (1.5) 



and summation over repeated (Greek) indices is implied. The first term in ( |1.2| ) describes free 
left- and right-moving electrons, whereas the second and the third terms describe forward 
and backward electron-electron scattering, respectively. The couplings g (> 0) and A (> 0) 
depend on the microscopic parameters of the lattice model and vp is the Fermi velocity. 
Normal ordering :: is carried out w.r.t. the filled Dirac sea. 



We shall treat the model in (IT) using the newly developed conformal field theory 
approach to quantum impurity problems by Affleck and Ludwig [|18| -|20|j. The basic idea of 
this method is to replace the impurity by a boundary condition, in the spirit of Nozieres' 
local Fermi-liquid theory of the ordinary Kondo effect 0. At the low-temperature fixed 
point, the long- wavelength properties are described by a conformally invariant boundary 
condition and boundary conformal field theory [[2l"H23| essentially determines all universal 
properties. Specifically, the theory predicts all boundary scaling operators that govern the 
asymptotic auto-correlation functions in the neighborhood of the impurity. As was realized 
by Nozieres many years ago || , the impurity response to an external bulk field is governed by 
the leading irrelevant boundary operators. Thus, knowing these, one can directly deduce the 
impurity critical behavior. The difficulty, though, is to identify the right boundary condition, 
although frequently the symmetries of the problem cut down the list of candidates to a small 
number. In short, each boundary condition is associated with a selection rule for combining 
the various degrees of freedom (such as charge and spin) at the boundary. The problem thus 
reduces to identifying the right selection rule. This can be done, according to the fusion-rule 



hypothesis of Affleck and Ludwig [|19|,f24]] by applying conformal field theory fusion rules ||25| . 

With only forward electron- impurity scattering present in (|1.3|) , the model can easily be 
cast on a form where boundary conformal field theory applies. Writing the Hamiltonian in 
terms of charge and spin currents (Sugawara construction [P5fl), the effect of the impurity 
is traded for a new boundary condition in the spin sector. At this point, however, two new 
elements enter the problem (as compared to the treatment of the single-impurity Kondo effect 
for free electrons). First, charge and spin excitations — although dynamically decoupled - 
are still connected via a selection rule of the type mentioned above. This must be carefully 
analyzed in the basis of states which diagonalize the interacting Hamiltonian. Secondly, by 
having both left- and right-moving electrons coupled to the impurity, left- and right-moving 
spin excitations are no longer separately conserved, only the total spin remains conserved. In 
technical terms (to be made precise), the chiral spin SU(2) x SU(2) symmetry of the critical 
bulk theory is not recovered at the fixed-point value of the electron-impurity coupling A. This 
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symmetry breaking introduces boundary operators with non-integer scaling dimensions, by 
a mechanism similar to that operating in the two-impurity Kondo problem |24| . 



The picture that emerges is consistent with that recently suggested by Furusaki and 



Nagaosa |13[ in their "poor-man scaling" analysis of the problem: the system renormalizes 
onto a low-temperature fixed point with the same critical exponents for impurity specific 
heat and susceptibility as in the two-channel Kondo problem |Tjl. Any asymmetry in the 



left-right electron- impurity coupling destabilizes this critical point, driving the system to a 
one-channel (Fermi-liquid like) fixed point. An interesting feature of our solution is that the 
leading-correction-to-scaling boundary operator at the symmetric fixed point is unique, in 
contrast to previous treatments of two-channel Kondo physics [|19j . Specifically, for vanishing 
electron-electron interaction, this implies a unique Wilson ratio also at low temperatures. 

To make contact with possible future experiments, one necessarily has to add electron 
back scattering off the impurity. This is so, since a "real" (Kondo type) spin exchange 

H K = X :*t (0)1^^(0) :• 5, (1.6) 

with ^ a (x) the electron field in (|1.4j), decomposes into 

Hk = T~Lf + T~~£b , (1-7) 
Tip being the forward scattering term in ( |1.3|) , and 

Hb = A [V1 )CT (0)|^^(0) ■ S + V^(0)K^(0) ' S ] ■ ( L8 ) 

The back scattering term Tig mixes left and right electron fields, and thus breaks both 
chiral spin (£77(2) x SU(2)) and charge {U(l) x U(l)) symmetry of the bulk critical theory. 
This results in the appearance of a relevant boundary operator which takes the system to a 
new fixed point, describing Kondo scattering in a Luttinger liquid. Turning off the electron- 
electron interaction, the system may still be represented by a two-channel Hamiltonian, but 
now coupled to a magnetic impurity in only one of the channels. This is known to give a fixed 
point theory with Fermi liquid exponents (as for the ordinary one-channel Kondo problem 
|3|j26|]). To include the effect of the electron-electron interaction is a more delicate problem: 
the mixing of left- and right-moving electrons in (|1.8|) obstructs a Sugawara construction in 
a basis where the interaction remains local. 

Not being able to attack the problem directly at a Hamiltonian level, we shall make 
the natural assumption that the full Kondo interaction may nonetheless be described by 
a renormalized boundary condition on the critical bulk theory. This is in accord with the 
expected behavior of any quantum impurity interaction, as discussed in [p0[| . Note that 



the relevant operator due to the Kondo interaction only couples to the boundary, i.e. the 
new fixed point in this scheme is a new boundary fixed point with the critical bulk theory 
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unchanged, and that all conformally invariant boundary conditions are scale invariant and 
correspond to such boundary fixed points. By demanding that the non-interacting limit is 
correctly reproduced, with analytic scaling in temperature for the impurity specific heat C imp 
and susceptibility Ximp, it turns out that conformal invariance together with the symmetry 
of the problem restricts the possible types of critical behavior to only two: Either the theory 
remains a local Fermi liquid in the presence of electron-electron interaction, or electron 
correlations drive the system to a new fixed point where (to leading order in temperature) 

C imp = Cl ((l/K p ) - l) 2 T(W)-i + C2 T and X i mp = c 3 T°, (1.9) 

K p = (1 + 2g/vF)~ 1 ^ 2 being the Luttinger liquid charge parameter, and 01,2,3 amplitudes 
depending on the scaling fields. The second case (|1.9|) agrees with the finding by Furusaki 



and Nagaosa referred to above |L3|]. In the conformal field theory scheme this scaling is 
implied by a novel selection rule for recombining the degrees of freedom at the impurity site, 
and we will discuss its properties in extenso. 

It should be noted that a faithful modeling of a magnetic impurity must allow for the 
possibility that the impurity carries a net charge giving rise to a screened local potential. 
In the case of a Fermi liquid its effect can be absorbed by passing to a new electronic basis 
with renormalized single-particle energies (which is the reason why potential scattering is 
often neglected in the ordinary Kondo problem). However, for a Luttinger liquid the effect 
of a local potential is more dramatic, as shown in ||. For this reason, potential scattering 
must here be treated on equal footing with the spin exchange interaction.^] Some attempts 
in this direction have recently been discussed in [Q, but we will not address the problem 
here. 

The rest of the paper is organized as follows. In the next section we introduce the 
Hubbard chain coupled to an impurity spin, and perform a continuum limit retaining only 
forward electron scattering off the impurity. In Section [ITI| we derive the finite-size energy 
spectrum and the corresponding spectrum of boundary operator dimensions, employing a 
particular variant of Affleck and Ludwig's fusion-rule hypothesis. This section also contains 
a matching of the Luttinger liquid selection rule for combining charge and spin excitations 



against that from a Bethe-Ansatz analysis of the Hubbard model P9|-pl|]. Employing the 
results for the boundary operator spectrum, the impurity critical behavior is identified in 
Section [IV| as that of the two-channel Kondo problem. In Section |V| we then consider 



the effect of adding electron back scattering off the impurity, thus treating the full Kondo 



1 Although a local potential is generated by the Kondo coupling A [27], it cannot simulate a realistic 
case — such as a magnetic ion trapped in a quantum wire — where one expects V S> A, V being 
the strength of the screened Coulomb potential. 
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interaction in a Luttinger liquid. Section [V]], finally, summarizes our results. Throughout 
the paper we try to provide sufficient information to make it essentially self-contained for a 
reader with some acquaintance with conformal theory. 



II. THE MODEL 

A. The Hubbard chain 

To make the physical picture clear we start by considering an explicit model for inter- 
acting electrons on a one-dimensional lattice, coupled to a single S = \ impurity. The 
Hamiltonian 

H = H H + H I (2.1) 

consists of a periodic Hubbard chain with nearest-neighbor hopping (t) and repulsive on-site 
interactions (U) 

= ~t J2^ C n + l,« + C Ll,*Cn,a) + U £ Vn,| > U > 0, (2.2) 
n n 

and couplings (J n ) of the electron spins to the impurity 

H I = Yl J n C lA a ^ C n,„- 5 > ( 2 - 3 ) 

n 

where we implicitly sum over repeated Greek indices as before. The electron creation and 
annihilation operators satisfy canonical anti-commutation relations 

and the number of electrons with spin a at site n is given by n nu = cl a c n a (without 
summation). The electron density n e is the expectation value of n n >. + n n , . 

Without interactions, ( |2.2| ) yields a free electron dispersion e(k) = — 2t cos ak, with a the 
lattice constant. The Fermi surface consists of the two points k = ±kp, kp = n e ir/2a. For 
small excitations (weak interaction) one may linearize the spectrum around the Fermi points, 
which gives rise to left- and right-moving particles with velocities =bup, Vp = 2at sin(7rn e /2). 
In terms of electron operators one has 

r — p -ik F na , ik F na In r\ 

C n,a — e C L,n,a + e C R,n,a > \ Z -°) 

with c Lna and c Rna referring to the excitations around k = —kp and k = kp, respectively. 
The sign of kp in ( |2.5| ) is a matter of convention; we use that c n<T is expanded in { e * fcna }. 
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These electron operators, defined on the lattice, are replaced in the continuum limit by the 
Dirac fields if) ra (x): 




2n 



(2.6) 



with normalization given by (|1.5| ). It is then straightforward to verify that the free part of 
the Hamiltonian (|2.2|) equals TixL i n (|l-2j ) with g = 0. 

As for the electron-electron interaction, it follows from substituting ( |2.5|) in Q2.2|) that 



U 



H 



t=0 



r,s=L,R 
C R,na C L,na C L,n-a C R,n, 



+ 



e eC R,n,1 C L,n,1 C R,n,l C L,n,l + ^.C. 



(2.7) 



describing forward, backward and Umklapp electron-electron scattering, respectively. Due 
to phase oscillations, however, the Umklapp term cancels away from half-filling (n e ^ 1), 
and in the continuum limit we recover the interaction part of (|1.2|) with g = Ua/2ir. In this 
paper we will consider gapless excitations only, which restricts us to U > and n e ^ 1. 

Finally, we rewrite the interaction with the impurity spin in ( |2.3|) . It interacts locally 
with a few sites and we may concentrate on n = and n = ±1, see Fig [I]. With J_i = J\ we 
recover in the continuum limit the previous electron-impurity interactions and ( |1.8| ) 

with the following couplings for forward and backward scattering against the impurity spin: 



Af = — (J + 2Ji), (2.8a) 
Ztt 

X B = — (J + 2Jicosvrn e ). (2.8b) 
2% 

Note that backward, but not forward, electron-impurity scattering depends on the filling 
factor when the impurity is coupled to more than one site.Q The first case of impurity inter- 
action addressed in this paper concerns forward scattering only, and from this construction 
it is clear that one can cancel backward scattering (As = 0) by coupling the impurity to 
the two nearest neighboring sites at quarter filling, i.e. J = and n e — \. For other filling 
fractions we may also fulfill A# = by allowing J 7^ 0. The next case of interest is the 
Kondo interaction, X F = \ B 7^ 0. It corresponds to coupling the impurity to one site only 
at arbitrary filling, i.e. J\ = in ( |2.8| ). 



2 There is a misprint in the corresponding formula given by Furusaki and Nagaosa jl3|| . Their 
expression for the coupling constant Jb is valid at half-filling only, whereas the rest of their analysis 
implicitly assumes other filling fractions. This does not affect their general results, though. 
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Let us mention that another way of canceling the back-scattering term is to couple two 
neighboring sites to S at half-filling and choose Jo = J\ (and all other J n = 0). With n e = 1, 
Umklapp processes now come into play in fl2.7j ), causing a mass gap in the charge sector 



The spin sector remains massless and describes a spin-- antiferromagnetic Heisenberg chain 
with two neighboring sites coupled antiferromagnetically to an impurity spin (Fig. |]a): 

7~(-spin — J Sn ' S'n+i + Jq{Sq + Si) ■ S , Sjy = Sq , (2-9) 

n 

with J = 4t 2 /U. This situation is similar to that considered by Eggert and Affleck [ Hfl , 
who studied a spin chain with two open ends coupled symmetrically to a single impurity 
(Fig. ^|b). On the basis of bosonization and numerical renormalization it was concluded that 
this system is in the same universality class as the two-channel Kondo model. The case above 
(Fig. |]a) where the spins coupled to the impurity interact mutually has been discussed by 
Clarke et al. . Also using bosonization, these authors proposed that two-channel Kondo 
behavior is manifest for this well. Their argument is based on a formal similarity 

between the bosonized versions of the Hamiltonians for the impurity-spin chain system and 
the two-channel Kondo effect, and the situation is somewhat less clear than for the open 
chain. It is therefore of interest to reconsider the problem, and we shall return to it below. 



B. Sugawara form 

In what follows we focus on the case away from half-filling (n e 7^ 1), described by the 
Tomonaga-Luttinger model Q1.2| ), together with the forward electron-impurity interaction 
( |1.3|) . It is convenient to rewrite the Hamiltonian in terms of charge and spin currents 

J r (x) = :^»^ r »: , (2.10a) 
J r (x)= :^{x)\<T a ^{x): , (2.10b) 

with r = L, R. These obey the (level-2) £7(1) and level-1 SU(2) affine Kac-Moody algebras 

H 

[jL(x),JL(y)] = ±4m6'(x-y) , (2.11a) 
[Ji(x), Jl{y)} = ie ahc Ji 2n5{x - y) ±—5 ab 5'{x - y) , (2.11b) 

R R R £ 

respectively, with J" the components of J r = (J^f, J-!, J*)- The normal-ordered products of 
the fields in ( |2.10| ) are defined by the usual point-splitting procedure: 

:WV(z)Vv,» = = Jim + 8)^ a {x) - + S)^ r>(7 (x))] . (2.12) 
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The terms of the Hamiltonian ( |1 . 1|) can now be identified with combinations of the quadratic 
forms : J r (x)J s (x) : and : J r (x) ■ J s (x) : . As an example, let us consider the case r = s = L. 
The normal ordering is again defined via point splitting, 



: Jl(x)J l (x) : = lim [J L (x + 8)J L {x) - (J L (x + 8)J L {x))\ , 

5— >0 

and we need to evaluate 

This can be done by using Wick's theorem and the Green's functions 

8 n 



with the result that 

J L (x + 8)J L {x) = :ipl >a (x + 5)ip La (x + 8)ijj\^{x)ijj L ^{x) : 



~i5 



(2.13) 



(2.14) 



(2.15) 



Hence, 



d 



:Jl(x)J l (x): = :^IA x ^lA x ^l,-A x ^l,-A x ) : + 2 '-AlA^^lA^ 



(2.16) 



(2.17) 



up to a total derivative from a partial integration. The analogous procedure for the spin 
currents yields (using that cr ail ■ cr vr} = 28 an 8 tw - 8 aii 8 vv ) 

3 3 c/ 

:J L (x)-J L (x): = -- ^[jx)tlJ L Jx)tljl_ a (x)tlJ Lt _ a (x)-.+- : ^ L>ir (x)i—^ L>ir (x) : , (2.18) 



and it is clear from ( 2.17 ) and Q2.18 ) that the first term of (|1.2| ) can be written entirely in 
terms of currents 



■AlA^^AlA^'- = \ ■■Jl{x)J l [x): +^ :J L (x) ■ J L (x): 



(2.19) 



(up to a total derivative). In the same manner, one straightforwardly replaces all other 
terms and arrives at the Sugawara form [25] of the Hamiltonian in (|1 . 1|) : 



Htl = ^ J dx | VF ^ 9 [: J L (x)J L (x) : + : J R {x)J R (x) :] 

+ F 9 [-Jl(x) ■ Jl(x): + :Jr(x) ■ J R (x) 



+ | [Jl{x)Jr{x] - 4J L (x) ■ J R (x)]} 



(2.20) 



and 
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Hf = \[Jl(0) + Jr(0)]-S, 



(2.21) 



with g = all /2tt for the Hubbard model. Rewriting the spin-current part of the Hamiltonian 
on matrix form : J r (x) ■ J r (x) ■ = \ ■ J r ,aix(x)J r ^ a (x) : where J r>a ^{x) = : ipl a {x)4> r ^(x) : 
— |^cr/x : i>$ ui^^rui^ ) '■ suc ^ that J r (x) = \cr afl J r ^^(x) , one notes that the spin currents 
are traceless (J r ^ a (x) = 0), and hence have no charge components. Thus, spin-charge 
separation is manifest in ( |2.20D . 

The spin interaction term J l(x) • J r{x) can be shown to be marginally irrelevant for 
g > [25] , and will be dropped henceforth.0 (The spin interaction is the only term in (|2.20|) 
that renormalizes, the remaining terms in H.tl being exactly marginal.) The piece of the 
Hamiltonian containing charge currents is diagonalized via the Bogoliubov transformation 

Jl(x) = cosh# j'l(x) - sinhfl jn(x). (2.22) 

The transformation is canonical, with the new currents jl also satisfying the U(l) Kac- 
Moody algebra 

[jL(x),j LR (y)] = ±Am5'(x-y). (2.23) 

Inserting ( |2.22| ) into ( |2.20| ), the transformation is found to diagonalize the charge Hamilto- 
nian when 

tanh2# = — - — . (2.24) 

v F + g 



Collecting the results, 



with 



n = n* TL + H F (2.25) 



+ | [:J L (x) ■ J L (x): + :J R (x) • J R (x):]} (2.26) 
the resulting critical bulk Hamiltonian, and with Tip defined in ( 2.21 ). Here 



v c = v FX 1 + — (2.27a) 
V vf 

v s = v F -g (2.27b) 



3 The term J l{ x ) • J r(x) is known to give logarithmic corrections to certain bulk correlation 
functions. Like any irrelevant bulk operator, however, it cannot influence the leading boundary 
critical behavior (see Sec. fTv| ). Hence, we can safely drop it. 
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with vf and g denned above. It is easy to see that H-tl ls invariant under independent global 
£7(1) and SU(2) transformations on the left- and right-moving (chiral) fields: 

VV^e^'W, t,^M ri/1 > r = L,R (2.28) 

with <p r a constant, and U r ^^ an element of SU(2). In fact, the Sugawara form of 'H* TL 
implies invariance under the larger chiral 

g = U(1) L x U(1) R x SU(2) ltL x SU(2) ltR (2.29) 

Kac- Moody algebra [33|] (with the subscript "1" denoting the level of the SU(2) algebras). 
However, the chiral spin symmetry gets broken by the impurity-electron interaction Tip for 
any value of the coupling A . This is in contrast to the Kondo effect for free electrons. In that 
case too there is only forward electron scattering off the impurity, but only with one type of 
chiral electrons (say, left movers). The critical bulk Kac- Moody symmetry is U(l) x SU 
(or U{\) x SU{2) 2 x SU(2) 2 in the case of two channels), and this symmetry is restored 
precisely at the strong coupling fixed point fll9 |. We return to this below. 

It should be emphasized that although the theory contains two distinct velocities, v c and 
v s , Lorentz invariance is manifest separately in the charge and spin sectors of 'H^l- This can 
be seen via bosonization [pS|. The terms containing charge currents represent a free boson 
theory, while the spin-current terms represent an SU(2) k = 1 Wess-Zumino-Witten model. 
Both theories are conformally invariant, implying also Lorentz invariance. This observation 
will be important for the applications to come. 



C. Boundary formulation 



At this point we reformulate the problem so that Cardy's boundary conformal field theory 
|2~TH2~3t applies. The boosted currents are defined in two-dimensional space-time (r, x), with 
the impurity sitting on the time axis. We can think of the time axis as a boundary with 
periodic boundary conditions imposed on the currents: 

j r {r, 0+) = >(r, 0_) , J r (r, 0+) = J r (r, 0_). (2.30) 

In Cardy's formalism no excitations (or two-momentum) may flow through the boundary, 
and hence periodic boundary conditions are excluded. To circumvent this restriction, we 
confine the system to the interval — i ^ x ^ i (taking £ — ► oo at the end), fold it in half, 

|34|| , see Figs. [3]a and |b. 



and x 



double the currents, and identify the two points x = 
The new currents, defined for x ^ only, are related to the old ones by 







= 3l(x) 


jI 




= 3r{~x 






= 3r(x) 


Jr 




= 3l{-x 



( J\ 



X 



JUx) 



J L (x) 
Jr(-x) 
Jr{x) 
Jl(-x) 



(2.31) 
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and the periodic boundary conditions in ( |2.30|) become 



ii(o) = m [ Ji(o) = j&o) 

j!(o) = ii(o) \ Ji(o) = Ji(o). 1 ' J 

Hence, a flow of excitations across x = in the original system corresponds to having them 
come in through one channel and then reflected back through the other. 

With this procedure, the critical bulk Hamiltonian is now defined on the positive x-axis 
only: 

n* TL = ^ E j/ x { J + '•&(*)&(*)'] 

+ | [: r L (x) ■ Jl(x) : + : J* R (x) ■ J R (x) :] } , (2.33) 

with the boundary condition Q2.32j ) imposed at x = 0. We also see that the electron- impurity 
interaction can be written as 

n F = X[Jl(0) + J 2 L (0)]-S. (2.34) 

The new Hamiltonian (^33|) is invariant under the full Q = U(l) x U(l) x SU(2) 1 x SU(2) 1 
Kac-Moody algebra and has conformal charge c = 4, i.e. we represent our c = 2 theory 
defined on — i ^ x ^ I by a c = 4 theory on ^ x ^ I. However, we may analytically 
continue the left-moving currents in (|2.33|) to the negative x-axisQ (Fig. Qc). From ( |2.32|) it 
then follows that these can be identified with the right-moving currents on the positive axis: 



ii(-z) = f R {x) f Ji(-x) = j 2 r (x) 



Hence we can formulate the theory in terms of left-moving currents only. This leads to the 
form of the Hamiltonian that we shall mostly use: 

H = H* TL + H F , (2.36) 

where 

H *tl = ^ E J £ dx {j ■■Mx)f L (x):+^ :Jl(x) ■ J* L (x)-] , (2.37) 



4 It follows from a Lagrangian description of the Dirac fields that the left-moving currents jl and 
J l are analytic in z = vr + ix, with v = v c and v s , respectively, whereas the right-moving currents 
are anti-analytic. 
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and with Tip given in (|2.34|) . The analytically continued currents satisfy periodic boundary 
conditions: 



j1M)=J'1W, Jl(-Q=JlW, i = l,2, (2.38) 

yielding a theory defined on a ring with circumference 11. 

Before proceeding, it is instructive to look at the special case of non-interacting (g = 0) 
electrons, i.e. with v c = v s in ( |2.37| ). For this case, an alternative construction is possible: 
Introducing two channels ( "flavors" ) of left-going fields, 

= ^L, a , rtjx) = -4 (2-39) 



the free part of the bulk Hamiltonian in (1.2) attains the form 



n = V ^j <lr:rU-r)i^A.M): (2.40) 

with the electron-impurity interaction written as 

H F = \ 1^(0)1^^(0):. 5. (2.41) 

With the simple transformation in (|2.39|) , one thus arrives at a chiral (left-going) represen- 
tation of the two-channel Kondo model ||. This has the Sugawara form fli"8| 

dx | - :J L (x)J L (x): +- :Jl(x) ■ J L (x): 

+ ^:J F L (x)-J F L (x):}+\J L (0)-S, (2.42) 

where the currents Jl (charge), J l (spin), and J F L (flavor) generate the affine U(l), S , f/(2) 2 
and SU(2) 2 algebras, respectively. 

The structure in ( [2.42j ) is not easy to obtain in the presence of electron-electron in- 
teraction. Technically, the construction requires all excitations to have the same velocity, 
which is the case only when g = 0. The result nonetheless suggests that the interacting 



problem exhibits two-channel Kondo behavior Since the forward electron scattering 

off the impurity only affects the spin sector, the fact that the electron-electron interaction 
pushes the spin and charge excitations apart by endowing them with different velocities 
seems irrelevant. However, to put the conclusion on firm ground, one must carefully check 
the role of the selection rule for combining charge and spin excitations when interactions are 
present. One should here recall that although the electron-impurity interaction is entirely 
in the spin sector, the charge sector may nonetheless contribute correction-to-scaling oper- 
ators, as in the (one-channel) Kondo effect for non-interacting electrons. A second reason 
for dealing with the electron-electron interaction "head on" is that it gives us an inroad 
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to attack the problem of Kondo interaction in a Luttinger liquid (by including backward 
electron scattering off the impurity (see Sec. |V|)). 

Returning to fl2.37| ), we see that the bulk Hamiltonian 7i^ L separately conserves the U(l) 
and SU(2) excitations in the two channels. This simply reflects the fact that the Kac- Moody 
symmetry is given by Q = U(l) x U(l) x SU (2)i x S77(2)i, as it must. The electron-impurity 
interaction, however, breaks the S77(2)i x SU(2)i symmetry. This appears similar to the 
effect of the Kondo interaction on free electrons. In a left- (or right-) moving description, the 
Kondo term breaks the single spin SU(2)i symmetry (or SU(2) 2 in the two-channel case) of 
the chiral electron Hamiltonian. However, at a special value of the Kondo coupling, XkoucIo = 
X* Kondo , the impurity spin can be "absorbed" in the electron spin current via a canonical 
transformation. The impurity-electron interaction disappears from the Hamiltonian, and the 
SU(2)i (or two-channel SU{2) 2 ) symmetry is restored (now generated by the spin current 
of the combined electron-impurity system). In this scheme, A^ ondo defines the local strong 
coupling fixed point. It is tempting to proceed in an analogous way for the present problem, 
and try to absorb the impurity via the transformations J l L (x) — > J l L (x) + S5(x) ,i = 1,2, 
judiciously choosing a special value of the impurity-electron coupling A*. However, these 
transformations are not canonical and couple the two spin currents at the impurity site: the 
SU(2)i Kac-Moody algebras for channel 1 and 2 are no longer independent. Thus, the full 
U(l) x U(l) X S77(2)i x SU{2)i symmetry of r H.* TL is not recovered at A*. Now, suppose we 
could rewrite the spin part of 'W rL in terms of the total electron spin current J(x) = J\ + J 2 L 
(in addition to some auxiliary degrees of freedom). The total current J(x) generates the 
diagonal subgroup of SU{2)i x SU(2)i, and, as may be easily verified, satisfies a level-2 
SU(2) Kac-Moody algebra 

[J a {x), J\y)] = ie abc J c {x) 2n5{x - y) + m8 ab 5'(x - y). (2.43) 

The impurity may now be absorbed without problem, using the single canonical transfor- 
mation J(x) — > J(x) + S5(x) = J'(x). The combined electron-impurity current J'(x) is 
conserved, and the SU(2)2 Kac-Moody symmetry is hence restored. 

To carry out this program, one needs a dictionary to translate from the U(l) X (7(1) X 
SU(2)i x SU(2)i formulation of the problem to another representation in terms of U(l) x 
U(l) x SU{2)2 x Q, where Q is the symmetry group of some auxiliary degrees of freedom. 
Fortunately, such a dictionary already exists in conformal field theory (the coset construction 
of Goddard, Kent, and Olive (GKO) |35|), and provides an elegant solution to the problem. 
Before exploiting it, however, we set up the formalism for studying the finite-size spectrum 
of the theory. 
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III. FINITE-SIZE SPECTRUM AND BOUNDARY OPERATORS 



The finite-size spectrum of a 1+1D scale- invariant theory provides important information 
about its critical behavior. This follows from a well-known result in conformal field theory 
|36|j : The energy levels in a finite geometry are directly connected to the (boundary) scaling 
dimensions of operators in the (semi-) infinite plane. More precisely, consider a conformally 
invariant theory defined on the strip {w = v r + ix \ — oo < r < oo, ^ x ^ £} in Fig. |]a, 
with r "imaginary time" and x the space coordinate. (The velocity of the excitations of the 
Hamiltonian is denoted by v.) Then impose a conformally invariant boundary condition, 
call it A, at the edges x = and x = £, and map the strip onto the semi-infinite plane 
{z = vt' + ix'\x' ^ 0} in Fig. [|b, using the conformal transformation z = exp(7nu/l) 
(implying boundary condition A at x' — 0). With E° the ground-state energy, one has 

E = E° + *-^, (3.1) 

where {E} is the spectrum of excited energy levels in ^ x ^ £, and {A} is the spectrum of 
boundary scaling dimensions in the semi-infinite plane. In a 1+1D quantum mechanical re- 
alization the boundary x' = coincides with the time axis (r' also being an imaginary time), 
and it follows that the boundary dimensions determine the asymptotic auto-correlation func- 
tions. In other words, for |r'| ^> x', 

(0(t', x')O(0, x')} - (0(r f , x f ))(O(0, x')) ~ (3.2) 

with O an operator with boundary dimension A. 

For the present problem two additional features appear, not present in the standard 
scenario discussed above. First, the chiral (here, "left-moving") Hamiltonian in ( j2.37| ) rep- 
resents a full 1+1D theory on the cylinder (via the folding procedure in (|2.31|) ): the second 



channel of left-moving currents simulates the presence of right-moving currents. Therefore, 
bulk dimensions appear in the finite-size scaling formula, disguised as sums of dimensions 
of left-moving operators labeled by the channel index. Secondly, 'Htl supports two kinds of 
excitations, charge and spin, with distinct velocities v c and v s (when g ^ 0). However, as we 
already noted, the charge and spin excitations are dynamically decoupled, and conformal 
invariance (including Lorentz invariance) holds separately in the two sectors. Summing up, 
one expects that (|3.1|) is replaced by 



E — E° = E c + E s — (E° c + E° s ) , (3.3) 

where 

E c - K = ^(Al + A*), E s - E° = ^(A] + A?) , (3.4) 
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{A^} being the boundary dimensions in channel j = 1,2 and sector a = c (charge), s 
(spin). (This structure of the spectrum has been exhibited in a Bethe Ansatz analysis of the 
Hubbard chain [2S|; cf. the following section.) 



As we have seen in Section |11 C[ it is convenient to represent a theory defined on a strip 
by a chiral theory on a cylinder. Formally, this follows from the vanishing of the energy- 
momentum tensor at the boundary |2TJ, implying that left- and right-moving operators 
coincide at the boundary. This is precisely what we used in (|2.35| ) when we continued the 
charge and spin currents to the negative x-axis. We know that the energy spectrum of 
this theory is in one-to-one correspondence to the boundary scaling dimensions due to the 
boundary conditions applied at the edges of the strip. However, these scaling dimensions 
are only a subset of all possible chiral dimensions of the bulk conformal field theory. One 
is thus faced with the task to pick out those chiral dimensions that represent the wanted 
boundary condition. Formally, this may be done by connecting the boundary condition 
to a selection rule that prescribes how charge and spin excitations are combined at that 
boundary. Knowing the selection rule, and allowing only boundary operators that preserve 
the symmetry of the Hamiltonian, the formalism unambiguously predicts the set of possible 
boundary dimensions. For a quantum impurity problem that we are dealing with here, 
the more intricate task is to identify the correct boundary condition (or, selection rule) 
that represents the presence of the impurity: According to a conjecture by Affleck and 



Ludwig ||r9| , |20| , any quantum impurity renormalizes, at the fixed point, into a particular 
conformally invariant boundary condition on the critical theory that carries the extended 
degrees of freedom (in our case, TCtl)- One way of identifying this boundary condition is to 
start with some known, trivial boundary condition on the critical theory, with no coupling 
to the impurity. The associated selection rule simply describes the allowed combinations of 
charge and spin excitations when there is no impurity present. Now, place a spin impurity 
at the boundary, and couple it to the electrons. By redefining the spin current as that of 
electrons and impurity, J(x) — > J(x) + S8(x) = J'(x), the electron- impurity interaction 
is removed at the fixed point. The spin quantum numbers {j} will be shifted accordingly: 
{j} — ■> {j'}. The new selection rule, describing the renormalized boundary condition, is then 
obtained from the old by substituting {j'} for {j}. The fusion-rule hypothesis by Affleck 



and Ludwig Jl9| suggests that the shift of quantum numbers is precisely governed by the 



conformal field theory fusion rules, in our case those of the S77(2) 2 Kac-Moody algebra: The 
set of states (conformal tower) labeled by a quantum number j is mapped onto new sets 
labeled by j', where 

f=\j-i\,\3-i\ + 1, -> min{j +1,2- j- §}. (3.5) 

This is the essence of the conformal field theory approach, to be exploited below. 

In this section we study the finite-size spectrum of Htli an d derive expressions for the 
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scaling dimensions {A^} on a form adapted to the impurity problem. We verify our result 
by matching it to that of the exact Bethe-Ansatz analysis of the Hubbard model. Bringing 
the electron-impurity interaction into play, we then use the coset construction to make 
a conformal embedding of the original SU(2)x x SU(2)i spin currents into SU(2) 2 x Z 2 . 
(This corresponds to writing the spin part of TI^l as a single SU(2)2 Sugawara Hamiltonian 
together with an Ising model.) With this proviso, we suggest a particular application of the 



fusion-rule hypothesis |19| , p3| , and absorb the impurity spin in the total electron spin current 
J(x), using the conformal field theory fusion rules for the SU{2) 2 Kac- Moody algebra. From 
this, the spectrum of boundary scaling dimensions in presence of the impurity spin is read 
off. 



A. Finite-size spectrum of TL 



TL 



Before applying these techniques to our Hamiltonian ( p. 371 ), we need to rewrite it in 
Fourier space. Introducing the Fourier-transformed currents 




1 

2^ 




dxe im i x JL ^\ (3.6) 



for i = 1,2 and m € Z, the mode expanded Hamiltonian takes the form 

n* TL =Y,K + K, (3.7) 



i=l,2 



where 



oo 

{t^o + ^E^-^ (3 - 8a) 



m=l 



1 9 1 

K = ^{g^O ■ Jh + 3 E J -m ■ J l ~ 24}- (3-8b) 

m=l 

This result is easy to verify. Consider first the U(l) currents, with 

:f L (x)f L (x) : = lim [f L (x + 5)f L (x) - (f L (x + 5)f L (x))] . (3.9) 

Noting that j\ is a dimension-1 analytic operator with argument z = v c t + ix, it follows 
from the operator product expansion |37j of f L with itself thatQ 



5 The factor of 2 in the singular term comes from the "level" of the U(l) currents, these being 
sums over a =1,1 fields. Cf. ( |2.16| ) for Jl(x) that satisfies the same algebra. 
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Iji(x + S)f L (x)) 



2 

~ 1)2 Lv 



where the last term is the correction due to the finite I. The inverse of (|3.6f) 



(3.10) 



- y 



m=—oo 



-im^x I Jin 



(3.11) 



together with (|3.9| ) and (|3.10|) , implies that 
1 



2tt 



c£e :f L (x)f L (x): = lim 



<5->0 



7T 



E 



2£ 



p "■" l " 7 1 7 1 -J 

^ JnJ—n 1 ro 

7T0 Z 



(3.12) 



Using the C/(l) Kac-Moody algebra ( 2.23| ) in Fourier space 



we may write 



j^ + 2^cos(n^)jl n j; + 2^ 



ne - in T 5 , 



(3.13) 



(3.14) 



n>0 



n>0 



so that the last term cancels the singular part of ( |3.12| ): 

2tt 



E 

n>0 



(3.15) 



Hence, ( p.8a|) follows from comparing (|2.37|) with (|3.12| ). The analogous treatment of 
the Fourier transformed spin currents, satisfying the SU(2)i Kac-Moody algebra 



r jia rfcbl 
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■ abc jic jrik , r 



ab rife 



(3.16) 



leads to ( |3.8b|) . As can be seen from the derivation, the Schwinger type terms in (|3.8a|) 
and ( |3.8b| ), including the constant ^, are due to the regularization of the spectrum (normal 
ordering of the quadratic currents). These terms encode the conformal anomaly [^] of the 
theory, and we shall discuss their role below. 

Given TC l c ( |3.8a| ) and TL\ ( |3.8b| ), it is now easy to extract the finite-size spectrum. Let 
us start with the charge Hamiltonian, 7i c = 7i l c + 7i 2 c , and make the connection to the 
original electron fields in the Tomonaga-Luttinger Hamiltonian (|1.2| ). By construction (cf. 
Eq. ( |2.31| ) and ( |2.35[ j) we can identify 



Jm JL,' 



Jm JR,'' 



(3.17) 



where 
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3^ m = iz I dxe^jdx) 



(3.18) 



are the Fourier transforms of the left /right-moving U(l) currents introduced in ( |2.22|) . Thus 



= cosh6J L)m + smh9J R>r 
j 2 m = cosh9 J RtTn + sinh6> J L>r 



(3.19) 



with 



j L =— dxe ±m ? x iipl (xU L (x): 



(3.20) 



the Fourier components of the original U(l) currents in ( 2.10a ). For later convenience we 
impose anti-periodic boundary conditions^ on ip ra (x), which then can be expanded as 



(3.21) 



This definition implies that the momenta for ip L are given by k = Tj( n + §)> so that the 

R 1 1 ~ 

"single-particle" energy levels ^v c k, w.r.t. the Fermi level k = 0, both satisfy 



e n = 5 ). 



(3.22) 



As follows from the definition of the integrand of (|3.20|) in terms of point splitting ( |2.12| ), 
we need the (finite-size) Green's function 



(±iS) £ y t 



(3.23) 



to obtain 



Jl = lim 



71 



2^-mO 
n+m ^ ^ 



(3.24) 



It is now convenient to introduce normal ordering in Fourier space. The filled Fermi sea 
occupies all levels for n < 0, and hence 



cr,r?i ' 



for r = L or i2. Using {VV )C r,m> ^,n} = ^r S <W5 m „ and 



(3.25) 



6 The condition Vv,o-( — £) = — ^V,o-(0 is consistent with periodic boundary conditions on the cur- 
rents and implies that the Fermi level lies between two energy levels in the finite-size spectrum. 
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n>0 



(3.26) 



we may finally write 



n+m ' 



(3.27) 



The interpretation of J^ >m is now straightforward: For m = it counts the net number of 
left(right) moving particles w. r. t. the Fermi sea, whereas for m ^ it excites particles m 
steps. 

Next we introduce a set of Kac-Moody primary states \ Ap) w.r.t. the charge and spin 
currents |57|| , defined by 



^,m|Ap) = 0, JL m |A P ) = 0, m > 0, 



(3.28) 



from which all other states can be generated. In an occupation number representation [38 

\Ap)=\Q Lh Q Ll ,Q Rh Q Rl ) (3.29) 

with the constraint 



Qh+ — Ql I [^ 1 . 



(3.30) 



Here | Q Ta ) denotes a non-excited state with Q TU the number of r = L, R electrons, carrying 
spin a, added to the filled Fermi sea. Hence, the total charge in channel r, Q r = Q r j + Q r i, 
is the eigenvalue of J r $ in (|3.27| ). Combining (|3.8a|) , (|3.19| ), and (|3.27| ), we thus have 



(3.3i; 



where we have introduced the eigenvalues of j'q, 



1 e e 
q2 =Q-±AQ— 



(3.32) 



labeled by the quantum numbers 



Q = Ql + Qr — Qli + Qli + Qri + Qri 
AQ = Qi — Q R = Q L <i + Ql^ — Q R ^ — Q R ^ . 



(3.33) 



As Ql, Q r G Z, it follows that Q, AQ G Z with the restriction that Q ± AQ is even. 

Turning to the spin Hamiltonian, H s , we proceed analogously. Writing J l - J\ = 3Jq z Jq 



in (|3.8b[) , and identifying Fourier modes, 



J, 



Iz 



J, 



L,rn 



j2z jz 

" R,m i 



(3.34) 
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the spin Hamiltonian ( |3.8b| ) and the Fourier transform of ( |2.10b|) imply that 



KIAp) = -JL{-(C^-Qiy--}\A P ). (3.35) 

(In obvious notation, Q 1 = Ql and Q 2 = Qr.) Introducing spin quantum numbers f = 
0, |, i = 1,2, connected to the particle numbers by 



4 



y w ■ -/ v-v'T ^J. 
the eigenvalues in (|3.35|) are expressed as 



f(f + l) = (Q\-Q\Y <1, (3-36) 



WH^jgW + l)-^}, J* = 0,i. (3.37) 
(This result can also be obtained directly by noting that J l ■ J l is the Casimir invariant of 



SU(2)i, with eigenvalues + 1) when acting on the primary states f33[. 



As a consequence of [H c , TC S ] = 0, the primary states factorize in charge and spin, 

\A P ) =|A P ) C ®|Ap) s , (3.38) 

and as the Hamiltonians are diagonal in channels 1 and 2, 

|A P ) S = Ij 1 )^ ®|j 2 ) (2) (3.39a) 
|Ap>e= ®k 2 ) (2) , (3.39b) 

the superscripts denoting the two channels. Note that the q l, s mix the original L and R 
channels and are defined by ( [3.32j ) in terms of the quantum numbers Q and AQ. 

To obtain a state of arbitrary particle number (Qlt, Qli, Qr\, Qr{) we undo the con- 
straints (|3.30|) by applying the operators J^-m — Jr-m ^ Jr-mt with m > 0, on | Ap). These 
operators flip spin within the left- and right-moving branches. The resulting states, call 
them | A), are still labeled by Qh a , which now may be any integer. Furthermore, by acting 
with the operators J r - m and J*_ m (m > 0) on | A), we also remove the constraint that | Q ra ) 
is non-excited, i.e. we allow non-filled levels below the highest occupied level: As is readily 
verified, the Fourier modes | J r ~ m ± Jr-m ( m > 0) create "particle-hole excitations" from 
the states | A) within each (5 m -branch. Any state of particle number (Qli,Qli,Qri,Qri) 
can thus be obtained from the set of primary states. (In the following, we label our states in 
the diagonal basis of TC, though.) As follows from ( |3.19| ) and ( |3.28| ), the primary states | Ap) 
are also primary w.r.t. to the diagonalized currents j l _ m and J % _ m and any other excited 
state is obtained by applying these operators for m > 0. 

To extract the spectrum including the energy levels of the descendant states just exposed, 
it is sufficient to note that the operators {j^ m ; m > 0, k = 1, 2} and { J k _ m ; m > 0, k = 1, 2} 
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act as "raising operators" w.r.t. the primary states in (|3.39a|) and (|3.39b|) . Explicitly, from 
(pi) and (13131) , 



(3.40) 



Thus, comparing with (|3.31|) , the descendant levels in (the diagonal basis of) the charge 
sector are obtained by adding energy to the primary levels in steps of ttv c / t. The resulting 
finite-size spectrum organizes into conformal towers of equally spaced energy levels, each 
tower having a primary level q l as base:[] 



TTV, 



+ K 

24 c 



Nl e N, i = 1, 2. 



with 



Q e -±/\Q e — 1 Q, AQ e Z. 



Similarly, ( |3.8b|) and ( |3.16D imply 



[Us , J k - m \ = ^rnJ 



(3.41) 



(3.42) 



(3.43) 



A comparison with ( |3.35| ) yields the finite-size spectrum in the spin sector, with two confor- 
mal towers (j l = 0, ~) per channel: 



itv. 



Nt 



Nl e N, i = l,2. 



(3.44) 



Putting Q = AQ = j 1 



and summing over the channels in (|3.41 ) and ( p.44 ), we 



obtain the ground-state energy E° = E® + E® with 



Ylt 



J = c,s 



(3.45) 



As is well-known, the finite-size correction to the ground-state energy E° = E\ + E\ of a 
conformally invariant Hamiltonian (defined on a ring of circumference L x ) scales as El = 

R 

—7tvc/12L x p9| , v being the velocity of the massless excitations. In the present case, we use 
a chiral formulation and must accordingly compare ( |3.45[) with E\ only. Putting L x = 21, 
we obtain c = 2 in the charge and spin sectors, respectively, yielding a total c = 4. This is 
what we expect, since we are using a c = 4 representation of our original c = 2 theory. 



7 The invariance of Tt l c implies that its eigenstates fall into the irreducible representations of affine 
U(l). The notion of a conformal tower is used interchangeably to denote these representations 
or the corresponding set of spectral levels. The analogous remark applies for 7i\, invariant under 
affine SU(2) 1 . 



23 



Given the finite-size spectrum, we now identify the corresponding scaling dimensions. 
Using (|3~4p , together with (|3.41| ) and (|3.45|) , we have 



K = \ (e'f " (-lye-^Y +K, i = l,2, NieN. (3.46) 



Similarly, ( |3.4|) , ( 3.44 ), and ( |3.45| ) imply for the scaling dimensions in the spin sector: 



A; = ^(/ + 1) + AC, f = 0,l 2 = 1,2, iV'eN. (3.47) 



The dimension of a composite operator is thus given by 

A = Al + Al + Al + Al (3.48) 



We recognize the scaling dimensions in (|3.46j ) as those of a U(l) Gaussian theory represented 



by free bosons with periodicity </> = </> + 2nTZ, with 1Z = e 9 / y2 pEOH - The same structure 
is hidden in (|3.47| ). Substituting Q ± = Q]- — ± (Q 2 — Q 2 ) for Q and AQ respectively, 



choosing e 9 = 1 and using (|3.36|) , (|3.46| ) gives (|3.47| ). For this special value of the periodicity 



(JZ = l/v^2), it is known that the symmetry is enhanced to 577(2) |4T| , and the scaling 
dimensions become those of an SU(2) k = 1 WZW model. 

It may here be worthwhile to recall the key elements of the analysis: Equation (|3.48| ) gives 
the spectrum of chiral dimensions of the analytically continued theory in the full complex 
plane. As we pointed out above, boundary scaling dimensions in the half-plane correspond 
to subsets of this spectrum, with one subset for each particular boundary condition. By the 
"folding procedure" in ( |2.31| ), the subset corresponding to the boundary condition ( |2.32| ) 



exactly coincides with the set of bulk dimensions of the original theory in the full plane. To 
pick out these dimensions from fl3.48| ) we must first identify the rule for assigning values 
to the quantum numbers (Q, AQ, j 1 , j 2 ) in presence of the boundary condition ( p.32j ). By 



the equivalence with the bulk problem, this is the same as identifying the selection rule 
governing the bulk energy spectrum of H^l- 

Assuming that the electron-electron interaction does not obstruct the choice of particle 
number (Qlt, Qli, Qnh Qui), (|3.33j ) and (|3.36|) imply the Luttinger liquid selection rule: 

j 1 = i(Q + AQ)modl 
j 2 = ±(Q-AQ)modl. 

Remember that Q ± AQ is even, which is consistent with the allowed values of j 1 . At this 
point we would also like to point out that this selection rule includes an implicit relation 
between the two charge channels 1 and 2: The definition of q l in terms of Q and AQ only 
permits certain combinations of the q l and q 2 conformal towers. In contrast, the conformal 
towers in the spin sector are not constrained by such a relation. 
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The scaling dimensions of the possible boundary operators are now obtained from (|3.48|) 



by feeding into (|3.46| ) and (|3.47|) the allowed values of (Q, AQ, j 1 , j 2 ) according to fl3.49|) . 



Not all of these operators appear, though, in the effective theory describing the scaling 
region. In general, this can be written as an expansion 

H = H* + Y,9kO k (0), (3.50) 

k 

where H* is the fixed point Hamiltonian, and gy. and Ok are the associated scaling fields 
and boundary operators. Corrections to scaling are produced by the irrelevant operators, 
and these must be invariant under the continuous symmetries of 7i (as must any relevant 
operators in absence of external perturbations). Applied to our case, TC = 7~(.TL + J2k 9k@k(ty, 
invariant under chiral U(l) x U(l) x SU(2)\ x SU(2)\. In other words, charge and spin are 
conserved separately in the two channels, implying Q = AQ = and j 1 = j 2 = for Ok- 
Hence, the only boundary operators appearing at the fixed point are descendants of the 
identity. This trivial content of correction-to-scaling operators could of course have been 
predicted directly from symmetry arguments, without invoking the finite-size spectrum. 



However, having the spectrum in hand, including the selection rule in ( 3.49 ), we can attack 
the more challenging problem of electron-impurity scattering. 

Before doing so, however, we shall compare our spectrum and selection rule to those from 
an exact Bethe Ansatz analysis of the Hubbard chain (of which r W rL is the fixed point theory). 
As was first shown by Woynarovich [^9| (see also p0| , pl| ), the Bethe Ansatz spectrum of a 
Hubbard model on a finite ring also organizes into conformal towers. Away from half-filling: 

E - E ~ ^(A+ + A-) + ^(A+ + A;) + 0{±). (3.51) 

Here N is the number of sites on the ring, v c and v s are model dependent charge and spin 
velocities, and Af and Af are given by 

= \ (± ± UD C + 5- )) 2 + Nt (3.52a) 



A ^ = i y Is ~i ±Ds ) (3 - 52b) 

with N^jNf" G N. As in the conformal approach, the positive integers label "particle- 
hole excitations" (although the notion of a "particle" or "hole" in a Bethe Ansatz basis is 
different from that used here). The parameter £ c is a non- universal function of the micro- 
scopic parameters, while (J c , I s , D c , D s ) are quantum numbers subject to the Bethe Ansatz 
selection rule 

= j(/„ + J.) (modi) 2Dci2ZW<eZ . ^ 
= ^Ic (mod 1), 
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To match our result in Q3.46| ) and ( |3.47|) to that in (|3.52|) , we note that the dressed charge 
£ c and our parameter e e both measure the strength of the Hubbard interaction. Explicitly, 
using the parameterization in Sec. |TJ, we have that e 9 = [1 + U/(2nt sin(a/ci?))] 1 / 4 , while 
£ c = \/2[l — U/ (87rt sin(afco))] [3D] , k playing the role of a "Fermi momentum" in the Bethe 
Ansatz formalism. Putting ko = hp, it follows that £ c = \[2e~ e + 0[(U/t) 2 }. To lowest 
order in U/t, A* + A 2 C and A] + A^ in ( 3.46| ) and fl3.47j ) indeed exhaust the same combined 
spectrum of scaling dimensions as A+ + A~ and A+ + A~ in ( |3.52 ). This can be seen from 
the following analysis. 

Let us first generate all combinations of scaling dimensions according to the Luttinger 
Liquid selection rule ( |3.49|) and show that they equal allowed combinations of scaling di- 
mensions according to the Bethe Ansatz selection rule ( |3.53|) . We may concentrate on 
N l c = iVJ = 0, as higher levels can be trivially mapped once the relation for = JVJ = has 
been established. Moreover, if we formally set Q = I c and use £ c = V2e~ 9 , we see that it is 
sufficient to check that AQ and 4D C + 2D S span the same range of integers when Aj + A^ 
= A+ + Aj. Whether this correspondence is possible or not depends on the selection rule. 

Consider Q = 2n, n G Z. By construction, AQ is even, and from ( |3.47| ) and fl3.49| ) it 
follows that for A\ + A^ = we have AQ = 4m \4m + 2\, m G Z, for n even \odd\. This is 
exactly reproduced in the Bethe Ansatz spectrum by choosing D s = (which is allowed as I c 
is even) and I s = n. According to the selection rule, the allowed values of D c are p \p+ |\, 
p G Z, for n even \odd\, i.e. 4D C + 2D S gives the same range of integers as AQ. The other 
possibility, A^ + A^ = |, is analogously combined with the above values of AQ shifted by 
2, which is reproduced by choosing D s — 0, I s — n + 1, and D c = p + | \p\, p G Z. 

Next, consider Q — 2n + 1, n G Z. Then A^ + A^ = | and AQ spans all odd integers. 
Choosing I s = n and D s = ±| is consistent with ( |3.53| ) and gives A+ + A^ = |. For I s 
even \odd\, the allowed values of D c are p + | p G Z. In both cases, 4_D C + 2D S spans 
all odd integers as required. This completes the first part of our comparison. 

In order to prove that the Luttinger liquid and Bethe Ansatz spectra are identical, we 
must also check that any other allowed combination of quantum numbers of the latter only 
reproduce energy levels already obtained. (Note that we only compare energy levels and 
not their degeneracies.) So far we have exhausted all combinations of I c and D c . By also 
allowing I s and D s to take any permissible value, one can show that all energy levels in the 
Bethe Ansatz solution for Nf- = N^ 1 = correspond to allowed combinations of Aj + A^ 
and Al + A^ with iV* = and iV* G N, respectively. The calculation is straightforward, 
so we leave it out for brevity. As before, it is trivial to extend the mapping to include 
N^ 1 , Ng G Z + , and thus the desired result follows. 



26 



B. Coset construction: Boundary operators in presence of Tip 



Having established the spectrum of Ti^Li we now include the electron-impurity term 

H F = \[J l L (0) + J 2 L (0)]-S (3.54) 

and explore how this interaction affects the boundary scaling dimensions. For this purpose, it 
is useful to rewrite the Hamiltonian 7i^ L +Ti.F in terms of the total spin current J = J l L + J 2 L : 
The chiral SU{2)\ x SU{2)i symmetry of r W TL gets broken by Hf, whereas the diagonal 
subgroup SU(2)2, generated by J, remains as a symmetry (cf. our discussion in Sec. [H|). 

In the conformal field theory formalism it is sufficient to work at the level of represen- 
tations of the SU(2)k algebras. Given a direct product of two irreducible representations 
(conformal towers) of SU (2)i, we thus ask how the states reappear using representations of 
SU (2)2- The answer is immediately obtained from the GKO coset construction [j35fl : Prod- 
ucts of two SU(2) 1 conformal towers (j x )i x (j 2 )\ decompose into products of SU{2) 2 and 
Ising model conformal towers (j)2 x (0) according to 

(0)i x (0)! = (0) 2 x (11) + (1) 2 x (e) 

(0) ! x (I)! = (i) 2 x (a) 

(1) ix(l)i = (0)»x(e) + (l) 2 x(l). (3.55) 

(Note that x (0)i is degenerate with (0)i x (|)i.) Here f — 0, \ and j = 0, |, 1 label 
the conformal towers of S77(2)i and SU(2) 2 , respectively, while (f> = II (identity), a (order 
parameter), and e (energy density) label the three conformal towers of the Ising model 
(c = |). The scaling dimensions in the SU (2)2 sector are given by 

A s = ^j(j + 1)+N s , N s eN (3.56) 

while 

( + N (1) 

Rising =1 ± + N (a) N G N (3.57) 
[l + N (e) 

are those of the Ising sector. 

All primary states in the new representation are now labeled by (Q,AQ,j,(j)), with the 
new selection rule obtained by combining (|3.49|) and ( |3.55|) : 



(0,11) or (l,e) Q even, \(Q + AQ) even 
(j.o) = { (0,e) or (1,1) Q even, \(Q + AQ) odd (3.58) 
^|,cr) Q odd. 
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At the Hamiltonian level, the conformal embedding SU(2)i x SU(2)i — > SU {2)2 x I sing 
implies that the spin part of 7Ytl decomposes into a sum of an SU{2) 2 Sugawara and a free 
Majorana Hamiltonian: 

rfx <^ - : J(x) • J(x) : +7] L (x)i—r] L (x) \ , (3.59) 
.1 < 4 dx J 

with 77.L a left-moving Majorana fermion. This enables us to absorb the electron-impurity 
term Tip = A«7(0) ■ S" into 7Y S by making the canonical transformation 

J(x) -> J(x) + Stf(a;) , (3.60) 

and choosing A = A* = In this scheme A* defines the (non-universal) fixed point coupling: 
the elimination of H.f by ( |3.60| ) restores translational (and thereby conformal) invariance, 
implying a fixed point theory. 

The redefinition of the total current in ( |3.60|) changes the rule for coupling conformal 



towers. Effectively, ( |3.60| ) adds an extra spin-| degree of freedom to the SU{2) 2 sector, 



leaving the U(l) and Ising sectors intact. As a result, the SU{2) 2 conformal towers (j) 2 get 
replaced, according to 

(0) 2 -(|) a , (|) 2 - (0) 8 or (l) a , (1) 2 -(§) a . (3.61) 



These are the SU{2) 2 fusion rules, see ( |3.5|) , describing the effect of combining a tower of 
spin j with a spin-| tower. Together with fl3.58|) this leads to the new selection rule 



I (0, a) or (1, a) Q odd, 

which governs the finite-size spectrum of r W TL +Hf, or equivalently, the spectrum of TH-tl 
with a modified boundary condition at x = representing the impurity spin. 

However, to obtain the scaling dimensions related to a boundary condition, we recall from 



the beginning of Section |TTT| that these are in one-to-one correspondence to the finite-size 
spectrum of the Hamiltonian with that boundary condition applied at both ends. In terms 
of the Hamiltonian ( |2.33 ), defined on x G [0,£], this corresponds to having impurity spins 



at both ends of the space interval. Passing over to the chiral formulation fl2.37p , we thus see 
that the boundary scaling dimensions of Ti, = 'Hp L + H,f are in one-to-one correspondence 
to the energy levels of the auxiliary Hamiltonian 

H' = H* TL + \[ dxJ{x)-[S 1 6(x) + S 2 S(x-£)}. (3.63) 
J-t 

The two impurity terms are removed at the fixed point by the canonical transformation 



28 



J(x) — > J(x) + S\5(x) + S 2 S(x 



(3.64) 



by which two extra spin-| degrees of freedom are added to the SU(2) 2 sector. The fusion 



hypothesis [IS] suggests that the effect is described by two repeated fusions of the spin 



j conformal towers with a, j = - tower. The new selection rule that emerges is readily 



extracted, using ( |3.61| ) twice to replace the SU(2) 2 towers in ( |3.58| ) by the resulting set of 
conformal towers. One finds, 



(0 or 1, II or e) 



Q even 
Q odd. 



(3.65) 



As before, AQ has the same parity as Q, i.e. they are both even or both odd. 

Summarizing, the spectrum of boundary dimensions in the presence of Tip is obtained 
from 



A = Al + A^ + A s + A 



I sing j 



(3.66) 



using (|3.65|) to insert the allowed combinations of (Q,AQ,j,(p) into the expressions for 



A* (i = 1,2), A s and A Ising in flOp , Qjjp , and (|3~57|) , respectively. 



IV. CRITICAL BEHAVIOR 

We now turn to the question of the impurity critical behavior, exploiting the result 
derived in the previous section. To do so requires two steps. First, the symmetry preserving 
boundary operators must be identified from the spectrum in (|3.66j ). With these in hand, one 
then selects the operator of lowest dimension to perturbatively calculate the leading finite- 
size corrections to the fixed point theory. Treating the (inverse) length as a temperature 
variable allows for the calculation of the finite-T scaling of physical response functions, such 
as the impurity contribution to the specific heat and magnetic susceptibility. 

It is tempting to make a shortcut and refer to the logic of our procedure to conclude 
that the impurity critical behavior must be that of the two-channel Kondo problem: The 
selection rule in ( |3.65| ), describing the renormalized boundary condition due to the impu- 



rity, descends from the Luttinger liquid selection rule in ( 3.49 ) via the coset construction 
and "double fusion". This selection rule in turn is nothing but the free-electron selection 
rule in disguise, "rotated" to a diagonal basis in the charge sector. One thus infers that 
the interacting problem cannot differ in essence, provided that the charge sector does not 
contribute a (non-trivial) leading correction-to-scaling boundary operator (which, by chiral 
charge conservation, is excluded). Although this line of reasoning is essentially correct, it is 
instructive to carry out the analysis "by hand" , and observe how two-channel Kondo physics 
emerges in the scheme proposed here. In fact, certain novel features appear: the leading 
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correction-to-scaling boundary operator (LCBO) is found to be unique at the fixed point, 
suggesting a universal Wilson ratio in the low-temperature, strong coupling limit. Also, 
by working out the solution explicitly in the new scheme we gain a number of important 
insights, to be exploited when extending the analysis to Kondo scattering in a Luttinger 
liquid (Sec. 0). 



A. Finite-temperature scaling 

To provide some background, let us briefly review the finite-temperature scaling approach 
to the problem, including the rather exotic mechanism that brings the irrelevant boundary 
operators to center stage. For more details, we refer to Refs. p^9| , |20| . 

Consider the (folded) system confined to a spatial interval [0, £] and in the presence 
of an external magnetic field h. By treating temperature T as an inverse length we 
may extract the low-temperature thermodynamics via finite-size scaling on a cylinder with 
circumference v(3. A convenient choice is to map the upper half-plane C + = {Imz > 0} 
(T = geometry in Fig. ^b) onto a cylinder T + = {w = (v/3/tt) arctanz} (finite-T geometry 
in Fig. [5|). At finite temperature the free energy separates into two pieces, describing the 
bulk and the impurity, respectively: 

F(p, h, A/) = 2£f bulk ((3, h) + f imp (f3, h, A/) , (4.1) 

with Xi the leading irrelevant boundary scaling field.0 By writing 2£ instead of £ in Q4.1|) , 
we let fbuik refer to the free-energy density of the original (unfolded) system; cf. Figs, ^a 
and |3]b and the remark after (|2.34 ) . 



From the standard finite-size scaling hypothesis f42fl , it follows that the reduced free 



energy density of a critical theory, defined on a cylinder with circumference v(3, scales as 

(3F/£ = f v(3 + Y{h(3)/vf3 + ..., (4.2) 

where the universal amplitude Yb u ik(0) = —nc/6 [^j. In our case, the critical theory is a 
sum of two conformal theories, each with c = 2, and with effective velocities v — v c and 
v s , respectively. Moreover, as the magnetic field h only couples to the spin sector, we may 
simply add the contributions from the two sectors: 



8 Note that we omit the irrelevant bulk fields in (|4.1| ). This can safely be done: If there is a 
correlation between irrelevant bulk and boundary operators, the corresponding term in a pertur- 
bative expansion of fimp in T can be shown to produce a subleading contribution compared to that 
coming from the leading irrelevant boundary operator alone. On the other hand, a bulk irrelevant 
field may produce nonanalytic corrections to fbulk, but this is of no concern to us here. 



30 



f Mk (J3, h) = E bulk - + ^—Y s {h(3) + . . . , (4.3) 

with E bulk a non- universal quantity and Y s (0) = —tt/3. Putting (3 = 1/T one thus obtains 
for the bulk specific heat: 

C bulk (T) = -T^f* = |(i + i)T + . . . , (4.4) 

with ". . ." denoting subleading terms as T — > 0. The bulk magnetic susceptibility can 
similarly be obtained from by expanding Y s to second order in h/T, and one finds 

Xbuik(T) = - d ^ ulh = — !— + . . . (4.5) 

In the same manner as for f bu ik one can write down a finite-size scaling ansatz for the 
impurity part of the free energy: 

f imp (P, h, A/) = £ imp + ^Y imp (hp, \iP y ) + . . . , (4.6) 

with E imp non-universal. The exponent y is an RG eigenvalue, connected to the dimension 
A of the leading irrelevant boundary operator by y — 1 — A. As the boundary scaling field 
Xj may couple to both the charge and the spin sector (A = A c + A s ), Y imp is not just a sum 
of a charge and a spin part. This will not be important, though, as we will later calculate 
fimp without any assumptions about Y imp . However, for Xj = there is only one scaling 
field h and it couples only to the spin sector. The /i-dependence of Y imp is therefore given 
by the same universal function as in a theory with a single velocity. 

From the scaling form of Yi mp one expects non-analytic temperature terms in fi mp for 
non- integer y and Xj ^ 0. Remarkably, the impurity specific heat and susceptibility vanish 
identically when A/ = 0, so these non-analytic terms become dominant. The mechanism 
behind this unusual behavior is particularly transparent within the conformal field theory 
approach: Since, at the fixed point, the impurity spin has been absorbed in the spin cur- 
rent J, the fixed point theory has lost all memory of the impurity. Specifically, the total 
magnetization 

ix{r(x)) = W(ftM, = 0) (47) 



0. (4.8) 



is insensitive to the impurity, and thus, from 

dY imp (h/3, Xj = 0) 
d(h(5) 

This implies that exactly at the fixed point 



31 



d 2 f 



dh 2 

and analogously 



d 2 Y imp (h(3, A/ = 0) 



WA, = 0) = -^-^ =-/3 - m Zr a :, ' =0, (4.9) 



/i=o d(h(3) 



h=0 



C %mp {h = 0) = -T^^ = 0. (4.10) 



To "put back" the effect of the impurity, the leading irrelevant boundary operator must be 
added, and this operator then produces the dominant term in the scaling in temperature. 

B. Boundary operators 

Following the procedure outlined in the beginning of this section, we first pick out the 
scaling dimensions in ( p.66|) corresponding to the symmetry preserving operators. 

It is easiest to start in the charge sector. Chiral U(l) invariance implies charge conser- 
vation in each channel, hence Q = AQ = 0. The selection rule in (|3.65| ) then implies that 
possible boundary operators can appear only in the products of conformal towers 

(Q,AQ,j,<j>) = (0,0,0,1), (0,0,1,11), (0,0, 0,e), and (0, 0, 1, e). (4.11) 



The first and fourth structures in ([4.11 ) contain only operators with integer scaling 



dimensions. As we shall see in the next section, these produce Fermi-liquid like, analytic 
scaling in temperature, and will be disregarded at this point. Turning to the second and 
third structures, we note that both contribute a relevant boundary operator, with A < 1.0 
The second possibility, (0,0, 1, 1), contains the spin-1 primary field, call it <fi, with A = \. 
However, the total spin is conserved, requiring all operators in the SU (2)2 sector to transform 
as singlets. This expels 4>, leaving only descendant operators in this sector with A ^ |. 
The third structure, (0, 0,0, e), contains the Ising energy density as a primary field, also 
with dimension A = |. Is it also suppressed by symmetry, or do we have to "fine tune" the 
theory to stay at the fixed point? Consider the latter alternative. The operator e originates 
from the conformal embedding of S77(2)i x SU(2)i into SU(2) 2 x Z 2 . The only parameter 
in Tt = Ti-TL + multiplying the SU (2)i currents is the spin velocity v s = Vp — g. A "fine- 
tuning" scenario thus implies that we can stay at the fixed point only for some privileged 
value of g, which, considering the known solution for non-interacting electrons, must be 
g = 0. Although not excluded a priori, this is not a likely situation. Let us instead explore 
whether there is any symmetry that suppresses e. 



9 An operator of scaling dimension A is relevant if D — A > 0. At a boundary, D = 1, hence 
A < 1 implies relevance. 
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The important symmetry to consider is invariance under channel exchange: 8 : 1 <-> 2, 
under which J\{x) <-> J 2 L (x). To test for 5, we connect e to the spin currents via the 
following observation: The difference of S77(2)i currents, J L (x) — J\(x), transforms in the 
adjoint (spin-1) representation of global SU{2) (with the Schwinger term subtracted): 

[J\\x) + J? (x) , Jl b (y) - Jl b (y)} = ie^J^x) - J 2 L c (x))6(x - y). (4.12) 

Knowing that the currents carry dimension 1, one is lead to the identification 

J\{x) - J 2 L (x) ~ <f>{x) x e(x). (4.13) 

This gives the correct assignment of spin (=1) and scaling dimension (A<^ + A e = § + | = 1). 
Now consider a channel exchange S. According to (|4.13 ), 



S : 4> x e -> -0 x e. (4.14) 

A consistent conformal representation hence requires that we assign evenness (oddness) to 
e {(f)) under £, or vice versa. Consider the first possibility: e even and odd. Since Tt is 
invariant under £, 4> and all its descendants become suppressed. The operator e on the other 
hand is allowed, and the required fine-tuning leaves us with a critical theory only at g = 0, 
according to the argument above. The first descendant of e, L_ie with scaling dimension 
I, here becomes the leading irrelevant operator. But this is exactly the boundary operator 



that drives two- impurity Kondo behavior [24 , in contradiction to the known two-channel 



behavior at g — 0. Agreement with established results in the non-interacting limit forces 
us to make the alternative assignment under £: e odd and <f> even. Now the e conformal 
tower gets suppressed, with no fine tuning necessary. Our analysis shows that there is no 
relevant boundary operator present: The only candidates appearing in (|4.11| ), (f> (A = |) 



and e (A = |) are expelled by symmetry. This yields the important conclusion that the 
theory flows onto a stable fixed point in the absence of external perturbations. 

The leading irrelevant boundary operator, call it Oi, is obtained by contracting <f> with 
the first SU(2) 2 raising operator J \. Oj = J \ • <fi, with dimension |. As expected, 
Oi is the same boundary operator that appears in the analysis of the two-channel Kondo 
effect, using the "charge-spin-flavor" scheme ||19|1 , reviewed in Sec. |T[ In that formulation, 



however, the flavor sector contributes a second operator with dimension |, call it Ofi avor . 
This is different from our approach where the leading correction-to-scaling operator Oi is 
unique. Consistency between the two schemes requires that the scaling field conjugate to 
Ofiavor vanishes at the fixed point. Based on the related case of "underscreening" , it has 
been argued that this indeed happens: When the impurity carries unit spin both the charge 
and flavor sectors contribute correction-to-scaling operators with the same dimension (=2) 
as the spin sector. Given a proper regularization of the theory, one can show that the 
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conjugate charge and flavor scaling fields — > at the trivial T > T% fixed point. There is 
evidence that an analogous mechanism is at play at the "overscreened" T < T% fixed point 
(with s = |), leaving only the operator Oj effective. Our result supports this picture, and 
allows for the definition of a universal Wilson ratio in the limit of vanishing electron-electron 
interaction. 



C. Low-temperature thermodynamics 



The low-temperature thermodynamics driven by Oi has been studied in ||19|| . In the 



next section we shall review this analysis, and extend it to theories with more than one 
velocity and boundary operators of arbitrary dimensions. Here we only quote the result for 
the impurity specific heat Ci mp and magnetic susceptibility x%mp- 



Ci m p{T, A/) 
Ximp(T, A/) 



A?9tt 2 



1 



Tln(— ) + O(r T), 



Xjl8 



ln( ^T ) + 0[(T ° T)C 



(4.15) 
(4.16) 



with r a short-time cutoff, playing the role of an inverse "Kondo temperature" . The scaling 
field A/ (conjugate to Oi) appears with the same power in fl4.15|) and (|4.16|) . Thus, given the 
known bulk response in (|4.4| ) and 
ratio 



, one may form the Xj- and T-independent Wilson 



Rw = 



Ximp/ X 



bulk 



bulk 



1 + 



1 - 



v F 



(4.17) 



In the g — > (v c , v s — > v F ) limit we recover the universal number | obtained by Affleck and 
Ludwig |19| for the two-channel Kondo effect. 

In the alternative charge-spin-flavor scheme for the two-channel Kondo effect, the flavor 
sector also contributes a dimension- 1 boundary operator. This produces a second term in 



Cimp, proportional to the square of a flavor scaling field A 



flavor^ 



To obtain the Wilson ratio 



R 



3 • 



one has to resort to the argument sketched above, suggesting that A 



flavor 



as 



one approaches the fixed point. In contrast, the reduction to one scaling field A/ is automatic 
in our approach. 

The results in ( 4.15 ) and ( [4.16 ) also apply to the associated spin chain problem, with two 
neighboring sites coupled antiferromagnetically with equal strength Jo to a spin-| impurity 
(cf. Fig. [|and the discussion in Sec. ||). This corroborates the result in [32], and shows that 
the closing of an open chain with its two ends coupled symmetrically to the same impurity 
does not affect the scaling behavior: both the open and closed chain exhibit the same "two- 
channel Kondo physics" ||14j| . This, from a naive point of view, is somewhat surprising: 
The two spins at the endpoints of the open chain are locked into an S = 1 state via the 
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interaction with the impurity. This causes overscreening of the impurity spin, resulting in 
the two-channel Kondo behavior. Upon closing the chain the endpoint spins get mutually 
coupled via an antiferromagnetic exchange, call it J. As J is increased above Jo one may 
have guessed that the local S = 1 state gets destabilized, thus inhibiting overscreening, 
and leading to a decoupling of the impurity spin ("Curie behavior"). This does not happen 
however. The results in ( J4.15| ) and ( 4.16| ) are valid independent of the relative strengths of 



the bare couplings J an d J: The low-temperature fixed point is stable, attracting a flow 
for arbitrary initial (bare) values of the couplings J and J . What does change, however, 
is the temperature scale at which one reaches this fixed point. For example, if the ratio 
Jo/ J is very small, the Kondo temperature will also be small, and one has to go to very low 
temperatures to see the "two-channel behavior" . 

D. Symmetry breaking perturbations 

It is instructive to study the stability of the forward scattering fixed point against sym- 
metry breaking perturbations. Consider first an asymmetric electron-impurity interaction, 

H' F = X L J L (0) ■ S + X R J R (0) ■ S , (4.18) 

with X L and X R arbitrary. This may be expressed as TL' F = Tip + 7~i per t where Tip is the 
symmetric electron-impurity interaction (with A = + X R )) studied above, and 

n pert = 1 -{\ L - \ R ) [Ji(o) - Ji(o)] • s. (4.i9) 



Tipert breaks invariance under channel exchange, and using ([4. 13|) one has 



H pert ~ [0(0) x e (0)] • S. (4.20) 

It follows that e enters as an allowed boundary operator, being contained in the spectrum 
( [4. 11| ). Having dimension A = |, it is relevant, and destabilizes the symmetric forward 



scattering fixed point .0 This is similar to adding a flavor anisotropic perturbation to the 



two-channel Kondo interaction, using a charge-spin-flavor scheme [26|. In this case the flavor 
sector contributes a dimension-^ operator which takes the system to a one-channel (Fermi 



10 Note that the spin-1 field cf) is odd under time reversal (T), and hence cannot appear alone, 
but only as a piece of a composite operator that respects T. As can be seen from ( |4.13| ), e is 
even under 7", and is thus allowed. Also note that exactly at the fixed point, the scaling field 
conjugate to e vanishes, and the invariance under e — > — e of the critical Ising model is recovered 
(Kramers- Wannier duality). 
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liquid like) fixed point. Again, equivalence between the two schemes requires e to perform 
the same way. An interesting question is how the scaling region about this one-channel fixed 
point is influenced by the electron-electron interaction. However, we shall not pursue this 
problem here. 

Instead, let us study the effect of adding electron backscattering on the impurity, that is, 
adding Hb in ( |1.8| ) to Hf- Hb breaks both chiral SU(2) and U(l) invariance. In particular, 
this means that the numbers of left and right moving particles are not conserved separately. 
As a consequence, AQ = Q L — Q R is no longer restricted to zero, and the charge sector 
can now make non-trivial contributions to the spectrum of boundary operators. The lowest- 
dimension operator with AQ ^ allowed by the selection rule in ( |3.65|) is obtained from 



(Q,AQ,j,0) = (O,±2,O,ll), (4.21) 

and has dimension A = \e~ 2e ^ |. Any \AQ\ > 2 will produce an irrelevant operator, and 
hence the composite primary operator from (|4.21|) is the unique relevant operator generated 



by breaking chiral U(l) invariance. (Note that the breaking of chiral SU(2) invariance is 
already effected by Hf-) As expected, Hb is thus a relevant perturbation that pulls the 
system away from the forward scattering fixed point, towards a new fixed point. The scaling 
behavior about this new fixed point is the topic of the next section. 



V. KONDO INTERACTION 

The problem we have studied so far (with only forward electron scattering off the impu- 
rity) is a model problem: we do not expect it to be observed in the laboratory. The reason is 
simply that by inserting a spinful (and chargeless) impurity into a, say, quantum wire, one 
necessarily produces a sharp scattering potential, the core of the interaction being a localized 
spin exchange. In contrast, the exclusion of electron back scattering off the impurity, i.e. the 
exclusion of large momentum transfers, requires an effective scattering potential that varies 
slowly in space. To have an experimentally relevant model we must therefore incorporate 
the back scattering term Hb, as in Eq. (|i~8f) , and study the full electron- impurity (Kondo) 
interaction TL F + 7i B in the strong-coupling, low-temperature regime. As we shall see, the 
insights gained from the treatment of the forward-scattering problem will turn out to be 
crucial in attacking this more difficult problem. 

Before plunging into the analysis, we briefly comment upon the recent work by Furusaki 
and Nagaosa [|l3j who exploit a version of "poor-man scaling" to analyze the problem. As for 
the ordinary Kondo problem |4c| , a set of scaling equations is derived perturbatively in the 



limit of weak electron-impurity coupling, and tentatively extended into the strong-coupling 
regime. This procedure suggests that the coupling increases indefinitely, implying a ground 
state where the impurity spin locks into a singlet with the conduction electrons, causing the 



36 



chain to break. A drawback of this method is that the scaling equations are formally valid 
only for weak couplings, and may miss out on possible intermediate-coupling fixed points, 
as in the two-channel Kondo effect. However, considering the analogy with the ordinary 
(one-channel) Kondo problem, it is quite likely that the result that ensues is in fact valid. 
Nonetheless, it is important to make an independent test of the result, using a method that 
is internally consistent. This is the aim of the following analysis. 



A. The non-interacting problem 

As a preliminary, let us look at the case of ID free electrons coupled to a localized spin 
(S — |) by the Kondo interaction Tix = Tip + H.B- The total Hamiltonian 7i = + Tix 
is given by 



H = y~ I <l-r 



dx 



(5.1) 



n 



K 



A 



K ■ 



A,M + 4,,(o)l k,,(o) + 1^(0) 



(5.2) 



Noting that only the even-parity part of the electron field couples to the impurity, it is 
convenient to pass to a Weyl basis via the canonical transformation 



^±A X ) 



V2 



^lA x ) ± ^rA~ x ) 



(5.3) 



From this construction it follows that both if) +gr and ijj- a are choral, left-moving fields with 
definite parity P under x — > —x: 

p^ ± A x ) = ±4>±A X )- ( 5 - 4 ) 

In the basis the total Hamiltonian takes the form 

Ho + ^ = ^ [dx ^tAx)i^ r A x )- +2A :< CT (0)K^(0): S. (5.5) 

71 r=+- •> X 

We recognize (|5.5| ) as identical to the Hamiltonian representing three-dimensional free elec- 
trons in two channels ("+" and "— "), coupled to a Kondo impurity in the "+"-channel 
only. This leads to a one-channel Kondo fixed point with a tt/2 phase shift of the single- 
electron levels in the "+" channel, the "— " channel being unaffected. (In fact, as pointed 
out already by Nozieres and Blandin (see also ||26|), the same conclusion holds for any 
channel-anisotropic Kondo interaction: the screening of the impurity is fully attained by the 
electrons in the more strongly coupled channel.) 
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The low-temperature impurity thermodynamics of the one-channel problem is that of a 
local Fermi liquid @. In particular, the impurity specific heat Ci mp and magnetic suscepti- 
bility Ximp scale as 

C imp ~T + 0(T 2 ), Xl m P ~T° + 0(T). (5.6) 

From (|5.5), this result also holds for one-dimensional free electrons coupled to a spin-^ 
Kondo impurity. 



B. The interacting problem 

The construction above is no longer useful when the electron-electron interaction is in- 
cluded, that is, when 




is added to T-Cq + T-Ck'- The interaction in (|5.7j ) mixes left- and right-moving fields, and hence 
becomes non-local in the Weyl basis. 

To make progress we must take a less direct route. We shall here exploit the expectation 
(see Sec. H and Ref. [pO|D that the full Kondo interaction Hk can be described as a renor- 
malized boundary condition on H^ L , analogous to the case of the "forward" interaction Hf 
studied in the previous sections. This is indeed a well-founded assumption: In the noninter- 
acting limit (g — 0) a canonical transformation on the even-parity spin current removes the 
impurity from the Hamiltonian in the Weyl basis (cf. Sec. |V A| ). This, "by construction" 
automatically leads to a change of boundary condition on the critical bulk theory. Turning 
on the bulk interaction, this boundary condition must still be present, although its effect 
(coded in the new selection rule for combining conformal towers) may change with a vari- 
ation of the bulk coupling. However — as we have just seen — when g ^ 0, Tix cannot 
be reformulated in terms of spin currents without violating locality of the electron-electron 
interaction. Therefore, we cannot identify the correct selection rule (boundary condition) 
by a redefinition of the spin current, as we did for Tip in Section |T|. In fact, we should not 
even expect that the new selection rule is simply related to the old one by a recombination 
of conformal towers in the spin sector only. To the contrary: Since the Kondo interaction 
H,k carries a charge component in the chiral basis, the charge sector is affected too. In other 
words, the selection rule for combining the two charge conformal towers may also change. 
This implies that the boundary operators may be composites of non-trivial operators from 
the spin and charge sectors. The fact that these are described by two distinct conformal 
field theories signals the novel aspect of the problem. 
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To obtain a sufficiently general framework for this new situation, we introduce a notation 
that does not make an implicit relation between the two charge towers, and denote a general 
combination of conformal towers by 



Here (Cj, D^) labels the U(l) channel-z tower (i = 1,2), while j and 0, as before, denote 
the SU(2)2 and Ising towers. The two U(l) towers are now treated as independent. In the 
description of 1~L* TL +Hf, as well as that of H.* TL with a trivial boundary, C\ = C 2 = Q and 
Di — D' 2 — AQ. These identities are the above-mentioned implicit relations between the 
charge towers and may be interpreted as part of the corresponding selection rule, whereas 
more general boundary conditions, associated with other types of impurity interactions, like 
TCr, may require C\ ^ C 2 and D\ ^ D 2 . Note that we still require Ci ± Di to be even for 
% = 1, 2, so as to preserve the spectrum within each conformal tower. 

The next step is to find the selection rule, and we first focus on the effect of Hk on the 
charge sector. As the quantum numbers Ci and Di can take any integer values, we assume 
that any selection rule for combining the U(l) conformal towers can be expressed by a linear 
relation: 



with a,/3, integers satisfying a5 — (3^ ^ 0. As we shall see, a few symmetry constraints 
severely limit the number of possibilities and leave us with only two permissible selection 
rules. Given these two rules for the U(l) x U(l) charge sector, we then consider all possible 
combinations of conformal towers from the SU(2) 2 x Ising sector. Each combination of 
states corresponds to a (composite) boundary operator, and those with scaling dimensions 
^ 1 are candidates for being the leading- correction-to scaling boundary operator (LCBO) 
that governs the critical behavior at the "Kondo fixed point" in a Luttinger liquid. As the 
Kondo interaction (|5.2| ) breaks the chiral U(l) (as well as the SU(2)) invariance of H.^ L , 
Ci and Di of the LCBO may now take non-zero values. However, global U(l) (as well as 
SU (2)) remains a symmetry, imposing other, weaker, constraints on these quantum numbers. 
(The conditions on the quantum numbers in the SU(2) 2 and Ising sectors are as before.) 
The list of candidate operators is then restricted by requiring that any LCBO respects the 
symmetry of the original Hamiltonian (including 7ix), and that the Fermi liquid scaling in 
( |5.6|) is correctly reproduced as g — > 0. (Note that a selection rule defines a boundary fixed 
point, and is valid for all values of the marginal bulk coupling g. Hence, given a selection 
rule, Fermi liquid scaling must emerge in the limit g — > 0.) 

With these preliminaries, let us make a first list of candidate LCBO scaling dimensions 



(C 1 ,D 1 ; C 2 ,D 2 ; j; 0). 



(5.8) 



Ci = a C 2 + (3D 2 + 77 , 
D x = 7 C 2 + 5D 2 +(, 



(5.9a) 
(5.9b) 
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at the TLtl + 7~Lk fixed point. From ( |3.46|) , ( |3.56|) , and ( |3.57| ), treating the two U(l) towers 
in (|3.46|) as independent, we have 



A = A c + A s , 



(5.10) 



where 



with 



A c = -{(gi) 2 + (g 2 ) 2 } +N C , iV c g N, 



qi = \Cie e + \D 1 e 
q 2 = \C 2 e e - \D 2 e 



(replacing the former q l of ( |3.32|) ) and 



A. 



°,1 








H 










n 




if 









+ N S , N s eN. 



(5.11) 



(5.12) 



(5.13) 



To cut down the list of possible scaling dimensions in ( |5.10| ), we next study the constraints 
imposed by the symmetries of the model. 



C. Symmetries and selection rules 

All states in the charge sector are combinations of states from the two U(l) conformal 
towers labeled by (C\, D{) and (C 2 , D 2 ), or q% and q 2 for short. Similarly, we label the Kac- 
Moody primary states of these conformal towers | q\) and | q 2 ), respectively. They transform 
under independent [/(l)-transformations as 

'/,)^^°' '/.,>• (5-14) 
It is convenient to introduce the linear combinations 



q = gi + q 2 = |(d + C 2 )e e + \{D X -D 2 )e 
Aq = qi -q 2 = \{C X -C 2 y + \(D X + D 2 )e 



and 



= | 
A4> = 



1 + 02 J 



2 Wl - 

so that we get a general combination of U(l) transformations as 



(5.15) 



(5.16) 
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I ft)® | g 2 ) =| gi ,g 2 > -> | gi, g 2 ) (5-17) 

with 

V = 9101 + 9202 = <?0 + AgA0. (5.18) 

In terms of these phase factors, global Z/ (1) is generated by </> and chiral U(l) by and A0. 
Hence, | qi,q 2 ) is global invariant if g = and chiral £7(1) invariant if g = Ag = 0. 
This is consistent with our previous notion of global and chiral U(l) invar iance in terms of 
Q and AQ, because q = Qe 9 and Ag = AQe~ d for the (bulk) Luttinger liquid and forward 
scattering problem, where Q = C\ = C 2 and AQ = D\ = D 2 . 

Consider now the original fixed point theory r W rL without impurity. Adding Tix breaks 
chiral, but not global, U(l) invariance. As 1~L* TL satisfies the Luttinger liquid selection rule, 
we see that the effect of adding the perturbation Tix is to break A<p invariance, i.e. to 
remove the constraint Ag = 0. Under renormalization the theory flows to a new fixed point, 
also associated with Tl^ L , but with a different selection rule, where, possibly, the two U(l) 
towers are decoupled. Whereas Q and AQ can then no longer be used to label the U{1) 
sector, q and Ag remain well-defined. Hence, at the new Kondo fixed point, the signature 
of adding backward scattering against the impurity is to allow operators with Ag ^ 0. 
The same arguments lead to the requirement q = 0, in order to preserve the global U(l) 
invariance of Tix- An LCBO must therefore not only be compliant with the selection rules 
( ^9|) , but also with 

C x = -C 2 and £>i = D 2 . (5.19) 

Next we consider the effect of a discrete symmetry. Although Tix breaks chiral U(l) 
symmetry, it is invariant under channel-exchange £ : 1 <->• 2, or equivalently L <-» R in ( |1.8| ). 
This translates into invariance under (g, Ag) — ■> (g, — Ag), which is equivalent to invariance 
under 

8: d «-> C 2 and D x <-> -D 2 . (5.20) 

Any candidate LCBO must respect this symmetry. In case g = Ag = 0, this is trivially 
fulfilled and A c = 0. We will return to this special case later and now focus on Ag ^ 0. 
Invariance under ( |5.2C| ) then implies that a LCBO must be a symmetric combination of 
operators with opposite signs of Ag. To force the coexistence of two such operators, we have 
to constrain the selection rules (|5.9| ) to be invariant under ( |5.20 ), i.e. they must satisfy 



C 2 = aC x -pDx + r], (5.21a) 
-D 2 = yC 1 - 5Dx + (. (5.21b) 
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If we now require at least one boundary operator, there must be a solution to fl5.9|), 
( g7TS| ) and ( gjTD - Inserting (pDD in QOal) and flOIap yields (1 + a)Ci = /3D 1 + 77 and 



(1 + a)Ci = /3I?i — 77, respectively, i.e. 77 = 0. Using (|5.9bp and ( |5.21b| ), we similarly get 
( = 0. We may therefore conclude that rj = ( = is a necessary condition for the selection 
rules. 

By demanding full consistency between (|5.9|) and (|5.21|) we may further reduce the list 



of possible selection rules: For instance, combining (|5.9a| ) and ( |5.21a|) requires (1 — a 2 )Cx = 
—af3Di + (3D 2 . Let us first consider the case a 2 7^ 1. Then a(3 = 0, as otherwise C\ would 
be a function of D\. (Remember that selection rules only give relations between conformal 
towers and should not pose constraints within.) (3 = implies C\ = 0, which is not allowed 
by the same reason. The other possibility, a = 0, implies C\ = (3D 2 . However, as C\ 
may be any integer, we can only allow (3 = ±1. Inserting this relation in ( |5.21b| ), yields 
+ 7)6*1 = 5D\. To avoid constraints within a conformal tower, we must then require 
(3' 1 + 7 = 5 = 0. Hence, the only solution for a 2 7^ 1 is a = 5 = and (3 = —7 = ±1. 
The next case, a = 1, implies (3Di = /3D 2 . If we assume (3 7^ 0, then D\ = D 2 , which gives 



C\ = C2 + (3Di using (|5.9a|) . But then (3 = 0, as C\ cannot be a function of Di, i.e. (3 ^ 
leads to a contradiction. The only possibility for a = 1 therefore is to require (3 = 0, which 
does not lead to a contradiction. The same result holds for a = — 1. 



Analogous treatment of ( |5.9b|) and ( |5.21b|) implies that either a = 5 = and (3 = —7 



±1 or else 5 = ±1 and 7 = 0. Hence, in total there are only six possible selection rules in 
the charge sector: 

d = C 2 and D x = D 2 (5.22a) 
d = C 2 and D x = -D 2 (5.22b) 
C x = -C 2 and D x = D 2 (5.22c) 
Cx = -C 2 and Dx = -D 2 (5.22d) 



and 



Cx = D 2 and D x = -C 2 (5.23a) 
Cx = -D 2 and D x = C 2 . (5.23b) 



The first selection rule (|5.22a ) is the Luttinger liquid selection rule that we start off with 
before we include Kondo scattering off the impurity. The effect of including Tix is to move 
us to a new fixed point, which may be described by any of the above six selection rules. In 
analogy to changing fixed point in the case of forward scattering off the impurity, or any 
other quantum impurity problem, we shall call such a transformation a fusion in the charge 
sector. It is a prescription for how the conformal towers are recombined when we change fixed 
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points. From ([5.22 ) and (|5.23|) it follows that there are six possible fusion rules that can be 
applied to the Luttinger liquid selection rule, and one of these should correspond to adding 
Tix to the fixed-point Hamiltonian. For instance, the fusion rule (C 2 ,D 2 ) — > (— C 2 , D 2 ) 
changes ( gggg) to (|5^D and (C 2 ,D 2 ) (D 2 ,-C 2 ) changes ( gggg ) to (^23^). Applying 
the correct fusion rule once to (|5.22a|) gives the new Kondo fixed point, and the selection rule 
can be used to extract the finite-size energy spectrum. Furthermore, we expect that applying 
the same fusion rule twice should give us the selection rule that determines the boundary 
scaling dimensions. It is easy to check that any fusion rule that takes ( |5.22a|) to any of the 
selection rules in ( [5 . 2 2| ) gives ( |5.22a| ) back after double fusion, whereas the fusion rules that 
take ( |5.22a|) to any of (|5.23|) give ( |5.22d|) after double fusion. We therefore conclude that the 
only possible selection rules for the boundary scaling dimensions are (|5.22a| ) and (|5.22d| ). 

We may now apply the symmetry constraint ( [5 . 1 9| ) to extract a "short list" of boundary 
scaling dimensions from the charge sector. The first rule, (|5.22a|) , together with ( |5.19| ), 
requires an LCBO to have C\ = C 2 = and D\ = D 2 an even integer. Using (|5.11| ) and 



A, 



We- 29 + N c 



P, N c e K 



(5.24) 



Similarly, the second selection rule, (|5.22d|) , combined with (|5.19| ), requires D\ = D 2 
and G\ = —C 2 an even integer, i.e. 



A c = ±pV* + N c 



P, N c e N. 



(5.25) 



The full boundary dimensions are obtained by coupling the SU(2) 2 and Ising conformal 
towers to the pairs of U(l) towers in ( |5.22a| ) and ( |5.22d| ), respectively. Starting with the 
SU(2) 2 sector, the j = \ tower is expelled by global SU(2) 2 invariance: Spin rotational 
invariance of the Hamiltonian 1~Ltl + *Hk implies that any LCBO must transform as a spin 
singlet, which, however, is missing from the j — 2 tower. Turning to the j = 1 tower, 
the primary operator <fi is expelled by the same reason. The lowest-dimension SU(2) 2 
singlet operator from this tower is J \ ■ (p. However, this is the same operator that drives 
critical scaling in the forward scattering problem. In particular, it produces a diverging 
impurity susceptibility as T — ► (cf. (|4.16| )), in conflict with the known Fermi liquid 
scaling (|5.6| ) in the g — > limit of the present problem.^ Assigning evenness to the j = 1 
conformal tower under channel exchange (as in the forward scattering problem) implies that 
J 1 ■ (f> is allowed by symmetry, and hence any selection rule must suppress this tower. The 



11 This conclusion still holds in the presence of non-trivial operator factors from the U(l) and/or 
Ising sectors. Whatever their dimensions, these can produce only T— independent constants in the 
scaling of the susceptibility with temperature (see Sec. |VE|) , 
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reverse assignment of parity under channel exchange instead implies that the j = 1 tower is 
suppressed by symmetry. Summarizing, the only possible contributions to an LCBO from 
the SU{2) 2 sector, consistent with established results for g = 0, are the identity operator 
and its descendants. We are thus left with the problem of gluing together the pairs of U(l) 
towers in ( |5.22a| ) and (|5.22d ) with the towers in the Ising sector. 

Let us start with ( |5.22aj ). Putting D% = D 2 = 0, the resulting identity towers can be 
combined only with the identity tower in the Ising sector (together with that of S77(2) 2 ). 
This is so, since the presence of a = cr or e tower would lead to a relevant boundary 
operator: the primary operators a and e both have dimensions < 1. However, at g — 
the fixed point is known to be stable (cf. Sec. [V A| ), excluding the presence of a relevant 
operator. In fact, this conclusion may be extended to g ^ 0: To remain at an unstable fixed 
point (that is, to maintain criticality) requires fine tuning of some parameter in the bare 
Hamiltonian. As g is the only tunable parameter in Ti.^ L (with a renormalized boundary 
condition replacing Hk), an unstable fixed point would imply non-criticality for all values of 
g 7^ 0. In other words, the total scaling dimension A = A c + A s of any boundary operator 
must be > 1. It is here important to stress that — by the raison d'etre of renormalization - 
any boundary operator allowed by symmetry will also appear at the fixed point. Therefore, 
one cannot argue that an LCBO with A > 1 can be obtained by forming descendants in the 
= a or e towers. If any of these towers were present, the corresponding primary operators 
of dimension A = ^ and |, respectively, would be present as well, implying an unstable 
fixed point. 

The lowest-dimension operators emerging from the identity towers are the first Kac- 
Moody descendants in the U(l) sectors, with A = 1: The marginal boundary operators 
1,2 (w) = j^ 2 (w) always appear in the charge sector, as the particle- hole symmetry of 
the original lattice model in (2.2) is broken away from half-filling.[^] Upholding particle- hole 
symmetry, the lowest dimension operators would instead be the second Virasoro descendants 
L_ 2 H in respective sectors, of dimensions A = 2. The next choice of U(l) quantum numbers, 
D\ — D 2 = 2 (p = 1 in (|5.24|) ), leads to a relevant boundary operator for any combination 
of Ising towers, and is therefore not allowed. In contrast, D\ = D 2 #5 4 (p ^ 2 in (|5.24j )) 
yields operators with A ^ 1 when combined with any Ising tower. Summarizing, the possible 



12 Although the Hubbard model ( |2.2| ) is not invariant under charge conjugation (particle-hole 
transformation) c n ,o- — > (— l) n e CT/ iCn iA1 off half-filling, its bosonized counterpart, the Tomonaga- 
Luttinger model ( |l.2[) , is invariant under the corresponding transformation ipr t a(%) — ► ^a/iipr^ix). 
The reason for this is that the normal ordering of the latter subtracts the symmetry-breaking 
terms. We may remove this accidental symmetry by explicitly introducing operators that break 
charge conjugation. 
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couplings of Ising conformal towers to the U(l) towers selected by ( |5.22a|) yield the following 
candidate LCBO dimensions: 

&LCBO = 1, \p 2 e- 2e + {0, i, |} , peN + 2. (5.26) 



Turning to the second selection rule for the Z7(l)xZ7(l) sector, (|5.22d| ), employing the same 



reasoning as above, one finds a second class of possible LCBO dimensions: 

^lcbo = 1, \e 2e + i, \p 2 e 2e + {0, i, 1} , pGN + 2. (5.27) 
Before exploring the critical behavior implied by the various dimensions in ( 5.26|) and 



( p.27p , two comments may be in order. First, note that no scaling dimensions of descendant 
operators — other than those of the identity — appear in ( |5.26| ) or ( |5.27 ). This is so, 



since an LCBO is the (composite) boundary operator with lowest dimension ^ 1, given 
a particular combination of conformal towers. Only if one, or several, of the non-trivial 
primary operators are expelled by symmetry can a descendant operator enter the stage (as 
in the forward scattering problem). This is not the case here. Secondly, the appearance of 
the e conformal tower requires the reverse assignment of parity under channel exchange, as 
compared to the forward scattering problem. That is, a consistent representation now forces 
the e tower to be invariant under channel exchange, while <f> and its descendants (already 
suppressed by the selection rules) change sign. (Note that there is no contradiction with the 
forward scattering the argument in Sec. [IV B| for the other assignment of parity no 

longer applies.) 



D. Impurity specific heat 

Analogous to the forward scattering problem, an effective scaling Hamiltonian 7i is ob- 
tained by adding a boundary term to the fixed point theory of Htl = + 'Hint- 

H = H* TL + A/0(O), (5.28) 

where O(0) is an LCBO with conjugate scaling field A/. By mapping the half-plane C + = 
{Imz > 0} for zero temperature onto the finite-T geometry (Fig. ||), T + = {w = v r + ix = 
(v/3/7r) arctan(z)} (with v = v c (v s ) for an LCBO from the charge (spin) sector^]), the 
partition function in zero magnetic field is written 



13 The finite-temperature geometry is defined in terms of r and x by T + = {— f3/2 ^ r ^ 0/2, x ^ 
0} with periodic boundary condition in r. When v c ^ v s , we represent it in complex notation by 
cylinders of different circumferences for charge and spin degrees of freedom, respectively. The only 
reason for introducing this complex notation is to simplify the derivation of correlation functions, 
whereas the calculation of the free energy is preferably done in terms of r and x. 
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so that Sf imp (/3, A/) = fimp(P, A/) - fi mp ((3,0) satisfies 



We have here used the decomposition in (|4.1|) and passed to a Lagrangian formalism, "tilde" 
and ( )t referring to T + . By a linked cluster expansion, 



=-jJ ^ dr(0(r, 0)) T " ^ // dr ^ d MO(n, 0)O fa, 0)} T , C + O(Af), (5.31) 

with ( ) TjC denoting a cumulant in T + . 

Here two cases must be distinguished: (i) O is Virasoro descendant of 11. Then (0(t, 0))t 
may be non-zero and hence the leading contribution to 5fi mp will be linear in A/, (ii) In any 
sector of the theory, O is Virasoro primary or a Virasoro descendant of an operator other 
than 1. Then (O(t,0))t = and Sf imp is quadratic in A/. To see how this comes about, 
consider a chiral (say, left-moving) Virasoro primary operator A(z) ^ II with dimension 
A (^ 0) in the half-plane C + . The expectation value of a chiral operator in a halfplane is 
the same as in the full plane, as translational invariance in one direction implies this in all 
directions. (Uniqueness of analytic functions implies that the expectation value is constant 
everywhere.) The scale transformation z — > z/a, where we later choose a = z, implies 

(A(z)) = z- A {A(1)), (5.32) 

and we conclude that (A(z)) = 0. (Note that the argument crucially depends on A being 
chiral: a non-chiral operator may pick up a non-zero expectation value in the presence of 
a boundary.) Any Virasoro descendant operator L_ ni L_ n2 . . . L- nM A{z) (rij > 0) has a 
vanishing expectation value as well, since (L_ ni . . . L- nM A(z)) = C- ni . . . £- nM (A(z)) = 
0, C-nj (j — l)---; n Af) being differential operators f44| . Now map C + onto T + . Again, 
by conformal invariance, (A(w))t = (^)~ A (A(z)) = 0. The transformation law for a 
Virasoro descendant is more complicated and relates one descendant in T + with a sum of 
operators from the same Virasoro tower in C + . However, as the expectation values of these 
are zero, it follows that the expectation value of a Virasoro descendant of A vanishes in 
T + as well. In contrast to this, descendants of the unit operator II may acquire non-zero 
expectation values in T + . As an example, the energy-momentum tensor is a descendant of 
II and satisfies (T(w))t = with c the conformal anomaly number and {z,w} = 

^~izj^m I ( ^dz/dw ) ^ ne Schwarzian derivative of the map C + — > T + . 

Let us first study case (i) where O is a descendant of 11. Since = dl/dz = 0, 

the energy momentum tensor T(z) = L_ 2 ll is the leading Virasoro descendant of 11 and 
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has A = 2. In the present problem, each sector contributes its own energy momentum 
tensor, T^z) = \ : j x L {z)jl{z) :, T 2 (z) = J : J 2 l (z)jI(z) :, T 3 (z) = \ : J(z) ■ J(z) : and 
T^z) = Tj S i ng (z), all with A = 2 and satisfying (Tj(z)) = in C + . In other words, there 
are four degenerate LCBOs for this case. Passing to T + , their contribution to the impurity 
specific heat is given by 

5f imp = - V -f / rfr(T,(r, 0)) r , (5.33) 



with 



(^H)t=^,W = |(^) 2 , (5-34) 

where c\ — — 1, C3 = | and C4 = | are the conformal anomaly numbers of the different 
sectors, and ui = v 2 = v c and v 3 = f 4 = f s are the corresponding velocities. Inserting (|5.34|) 
into (|5.33| ), integrating, and summing over the sectors, it follows that 

6 ^-f f • 

producing a linear specific heat 

j=i i 

where we have used ( 4.10p . This is the dominant contribution to Ci mp that is linear in the 
scaling fields: higher-order descendants of II produce higher powers in temperature. 

We now turn to the more interesting case (ii) where the candidate LCBOs are not Vira- 
soro descendants of the unit operator. The one-point function in (5.31) then vanishes, and 
we are left with 

A 2 r r p/2 

Sf imp = ~ d n dT 2 {0(n, 0)O(t 2 , 0)) T + O(Af). (5.37) 



2/? J J -13/2 

According to our symmetry analysis, any LCBO that is not a descendant of II must be a 
Virasoro primary operator.^ It is therefore sufficient to consider the case when O is Virasoro 
primary, for which the two-point function in C + takes the familiar form 

(0( Zl )0(z 2 )} = A (5.38) 
[Zi — z 2 ) 



14 All Kac-Moody primary operators allowed by symmetry, as well as the first Kac-Moody descen- 
dants in the charge sector, are Virasoro primary. 



47 



Here A 0) is the scaling dimension of O, and A is a normalization constant (to be 
determined). Using z — > w = v r + ix = (v/3/ir) arctanz and the transformation rule for a 
Virasoro primary operator, we can get (0(wi)0(w2))t- However, as O may be composed 
of operators from sectors with different velocities v, we must perform the transformation in 
each sector independently. The expression simplifies somewhat on the boundary x = and 
becomes 



(O(n,0)d(r 2 ,0)) T = — - A — . (5.39) 

|fsin(|(n -r 2 ))| 2A 



J c s 



where A = A c + A s , and the subscripts refer to the charge and spin sectors, respectively. 
The integrand ( 5.39| ) of ( |5.37| ) is even and periodic, with period /?, and can be replaced by 



a single integral over r = T\ — r 2 . Putting u = tan(7r//3), and inserting a short-time cutoff 
t = ef3/ir, yields the expression 

" = "w^F? ^ X»« >^ (5 ' 40) 

Let us first study the case when A ^ 1 is an integer. We can then make the expansion 
(1 + u 2 ) A ~ l = 1 + (A - l)w 2 + . . . + m 2 ^" 1 ), which yields 

Aane « 2A 2A - 1 (tane) 2A_1 \ 2A-3 1 ; K ' 

(5.41) 



In the limit of small e, i.e. r T — > 0, equations Q5.40|) and (|5.41|) give 

Cimr = -T^ = ( J B l )2 3 (2 f_3) ^°- 2Ar f 1 + °K*0D ' ( 5 ' 42 > 

The case of non-integer dimension A > 1 requires a more lengthy analysis. One partial 
integration of I yields 

' ~ ^t-^r (^Cfr.o^ri), ^ ,5,3, 

with Ii defined by 

00 , (l + u 2 )*-'- 1 



h = / du ^ > , j G N. (5.44) 
When A < |, /i is finite as e — > 0, while A = | produces a logarithmic singularity: 

/, - { * B( ? ~ A '^' A < f as r T - (5.45) 
In — -7= + constant, A = £ 
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with B(p, q) = T(p)T(q)/T(p + q) the Beta function. 

For A > |, I\ diverges algebraically with the cutoff, and we must perform a second 
partial integration to identify the rate of divergence. This leads to 



(5.46) 



with I2 defined as in (5.44). In the interval | < A < |, 1% is finite as e — > 0, while for A = | 
one gets the same logarithmic singularity as for A = |, i.e. ^(A = |) = ii(A = |). The 
procedure is now iterated, J J+ i being the remainder from partially integrating Ij. However, 
the terms in T thus generated are always subleading compared to those coming from the 
partial integrations of / and Ii in ( |5.43| ) and ( |5.45|) , respectively. Thus, collecting the results, 
we obtain for the impurity specific heat 

2A „2^3-2A T , A — 1 ™- A \ 3 



7r 2 r d-2A T + _ _ _ A = 1 or A > 



\2a 3(2A-3)" - ^ — ^ 2 

C *» = 73331 x 2 < A " i^^d " A - 5> T2A " + ■ ■ ■ 1 < A < I (5.47) 



where ". . ." denotes subleading corrections. 

Combining (|5.36|) and (|5.47|) with the result from the symmetry analysis, ( |5.26 ) and 
( p.27| ), we find that there are only two distinct possibilities for critical scaling: 

(i) C imp = 0(T), (5.48) 

and 

\2a 1/K p 1 

(ii) C imp = -^^(t - l) 2 B(l - 5^, |) T~p + O(T), (5.49) 

with K p = e -26 * = (1 + 2g /vf)~ 1 I 2 the Luttinger liquid charge parameter. The first case, 
( 5. 48] ), is implied when the LCBO carries dimension A = 1 or A > |. In contrast, the 



leading term in the second case, (|5.49j) , is driven by a composite LCBO of dimension A = 
A c + A s = \{e 2e + 1), corresponding to the second entry in ( |5.27| ).P| In terms of quantum 
numbers (Ci, D\\ C2,D 2 ; j; (ft), the operator is given by the sum of 

(2,0; -2,0; 0; e) and (-2,0; 2,0; 0; e). (5.50) 

Hence, the charge sector contributes a channel-symmetric combination of primary operators 
with q = and Aq = ±2e d (ci. discussion after ( |5.20|) ), which can be explicitly expressed 



in terms of vertex operators of free boson fields. The contribution to the scaling dimension 



15 Since we assume the interaction to be weak (g <C 1), we take \(e 26 + 1) < 3/2. 
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is A c = ^e 29 = l/2K p . Only by combining these operators with the Ising energy density e 
does one obtain a boundary operator of dimension A ^ 1, as required for an LCBO. The 
linear term in ( |5.49| ) comes from subleading terms generated by the same operator, as well 
as from leading terms due to the marginal operators j 1 (z) and j 2 (z). As we have seen, the 
latter operators are always present, in case (i) as well, due to the breaking of particle- hole 
symmetry. 

The specific heat in Q5.49 ) exhibits the same anomalous scaling in temperature as found 
by Furusaki and Nagaosa [|l3j. Also, the way the anomalous term vanishes as K p — > 1 (t? — > 0) 
is identical to that obtained in Although more work is needed to firmly establish which 

of the two cases applies, ( |5.48| ) (Fermi liquid) or (|5.49|) (non- Fermi liquid), the second, non- 
Fermi liquid case is clearly favored considering its emergence in an independent analysis. 
However, a caveat is advisable. In a recent study, Schiller and Ingersent treats a simplified 
model of a magnetic impurity in a "reduced" Luttinger liquid, composed of right-moving, 
spin-up electrons and left-moving spin-down electrons ||45|| . This problem, with only two 
branches of electrons (compared to the four branches of the full problem: two chiralities 
with two spin projections each) is mapped exactly onto the Kondo effect in a Fermi liquid, 
with Ci m p as in ( |5.48|) . Although the relevance of this simplified model to the Kondo effect 
in a full Luttinger liquid remains unclear to us, the result may be taken - as argued by the 



authors in ||5| - to give some indirect support to a Fermi liquid scenario. In any event, it 
is reassuring that our exact analysis gives room only to those two scenarios that have been 
conjectured in the literature. 

Before closing this section, let us point out that the impurity specific heat ( |4.15|) for the 
forward-scattering problem follows from the third case in Q5.47|) by inserting A c = 0, A s = |, 
and the value of the normalization constant for this case, A = (J-i ■ 4>\J-i • 4>) = 9.Q The 
amplitude in ( |5.49 ) can similarly be calculated by evaluating the norm of the corresponding 
LCBO. However, as we shall find in the next section, the favored case, with the LCBO in 
( |5.50| ), does not lead to a universal Wilson ratio when g ^ 0. For this reason we here leave 
the amplitudes undetermined. 



16 The charge parameter is defined as K p = [(l—g2/^VF)/(^+92/'^VF)] 1 ^ 2 in |l3|]. For small 52 > and 
switching from the conventional normalization of Fermi fields, {ifj r>a (x), ipl,fi(y)} = S rs 5 afJ- 5(x — y), 
to our normalization in ( |1.5[ ), this is the same as our definition of K p with g2 = irg. 

17 Note that (J_i • <p\J-i ■ (j>) = {(j) a \ JfJ-i I 4> h ) = (4> a I [ J x > J -i\ I 6 )> where the second identity 
follows from I (f> b ) being a Kac-Moody primary state. By exploiting the SU (2)2 Kac-Moody algebra, 
and the fact that </> transforms as a spin-1 object under global SU(2) (generated by Jo), one 
straightforwardly arrives at the result, A = 9. 
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E. Impurity susceptibility 



In the presence of a magnetic field, the partition function takes the form 

r/9/2 



e 



The shift of the magnetic susceptibility due to the impurity, 



d 2 f- 

U J VI 



imp 

Ximp — Q^ 2 

may thus be expanded to second order in Xj as 

"P/2 



h=0 



(5.52) 



Ximp = 4tt^ // dx ^ dx2 J-J dn. . .dT 3 (J z (T 1 ,x 1 )J z (T 2 ,x 2 )O(T 3 ,0)} TtC 

X 2 r r°° r r^/ 2 

+ 8^5 J J dXldX2 y---y^ /2 rfr i--- rfr 4(^( r i^i)^( r i,a ; 2)O(r3,0)a(r4,0)) TiC . (5.53) 

Let us consider the non-Fermi liquid scenario with O given by ( |5.50| ). The U(l), SU(2) 2 
and Ising sectors are decoupled, and hence there are no dynamical correlations between 
operators belonging to different sectors. As O contains only the identity as SU{2) 2 factor, 
it follows that 

( J z (r h x x )... J z (r m , x m )0(r, 0)) TjC = , (5.54) 

using that the LCBO is Virasoro primary, i.e. (0(w))t = 0. The expectation value of J z 
vanishes by the same reason, and we may decompose 

- (J*0 3 ) (J z 2 O a ) - (J z 4 ) (J z 3 ) = 0. (5.55) 

Hence, we infer that O does not give any contribution to Ximp to O(Xj). The same conclusion 
holds for any candidate LCBO obtained in Sec. [VQ , as all of them are Virasoro primary 
and contain only the identity as SU(2) 2 factor. Higher order terms in an expansion of 
Ximp in Xi can similarly be shown to vanish. The leading contribution to Ximp is instead 
given by the lowest- dimension boundary operator that has a non-trivial SU(2) 2 factor. By 
our symmetry analysis in Sec. |V G , this is given by the SU(2) 2 energy-momentum tensor 
T 3 (z) = | : J(z) ■ J{z) :, of dimension A = 2. Analogous to the ordinary Kondo problem 
T 3 (z) produces a finite impurity susceptibility (to first order in Xj), 

Xi mp ~T° + 0(T). (5.56) 



Thus, comparing with the free case ( |5.6|) , the electron-electron interaction is seen not to 
influence x%mp- the impurity remains completely screened, in agreement with the result of 
Furusaki and Nagaosa [13|. Also note that, by ( |5.49| ) and ( |5.56|) , the favored no n- Fermi 
liquid scenario implies a Wilson ratio that is non-universal, and depends on temperature. 
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VI. SUMMARY 



We have studied the low-temperature properties of a spin-| magnetic impurity coupled 
to a one-dimensional interacting electron system. By turning the problem into a bound- 
ary critical phenomenon and using conformal field theory we have reached the important 
conclusion that the symmetry of the problem admits only one of two possible fixed points 
describing the local electron-impurity composite: Either the theory remains a local Fermi 
liquid in the presence of electron-electron interaction (as for the ordinary Kondo problem 
with free electrons) or electron correlations drive the system to a new fixed point with an 



anomalous specific heat, identical to that proposed recently by Furusaki and Nagaosa [JT3 
We have also shown that the suppression of back scattering off the impurity destabilizes both 
fixed points and produces an impurity critical behavior identical to that of the two-channel 
Kondo model, but with a new Wilson ratio. 

The non- Fermi liquid fixed point is distinguished by the presence of a leading-correction- 
to-scaling operator of dimension A = A c + A s = |e 261 + |, corresponding to the energy level 
E', withQ 

E'-E = ^e 2e + -f\. (6.1) 

By numerically computing the finite-size energy spectrum of the two-impurity auxiliary 
problem and checking for the presence of E', one has, in principle, a diagnostic tool for 
deciding which scenario is realized: Fermi- or non-Fermi liquid. Unfortunately, with our 
approach we cannot derive a unique and complete finite-size spectrum at the non-Fermi 
liquid fixed point as the selection rule in ( |5.22d|) only applies to the charge sector; there are 
still a multitude of ways of coupling the two U(l) conformal towers to those in the SU{2) 2 
and Ising sectors. When deriving possible LCBO dimensions we were helped by symmetry 
constraints, which, however, are no longer applicable when considering the full spectrum. 
This fact makes the check against numerics more difficult, as one is essentially restricted to 
search for the single level E'. 

To make contact with future experiments clearly requires a more complete description of 
the system, as well as the inclusion of potential scattering off the impurity (see Sec. H). In the 
case of a magnetic defect implanted in a quantum wire, the observable of greatest interest 
is the shift of the average conductance due to the impurity, as this is the quantity most 
easily accessible in the laboratory. Other important characteristics include the local spin 



IN 



Note that we here refer to the spectrum of an auxiliary problem with two impurities, one at 



each end of the interval [0, £]. Cf. the discussion on double fusion in the text after Eq. ( 3.62 ) and 
(|5T23bD. 
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and charge Green's functions, the scattering matrix, and the residual entropy. Considering 
the success of conformal field theory techniques for obtaining these quantities in the multi- 
channel Kondo problem |19|| , we judge that the approach as presented in this paper will be 
equally powerful. We hope to return to these, and related issues in a future publication. 
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FIGURES 

FIG. 1. The Hubbard chain interacting with an impurity spin S at a few sites. Solid lines 
represent electron hopping and dashed lines spin interactions. 

FIG. 2. The Heisenberg antiferromagnet in ID coupled via two neighboring sites to an impurity 
S. In (a) the chain is closed, in (b) open, although linked via the impurity. 

FIG. 3. Equivalent representations of the Luttinger liquid. The original c = 2 system (a) 
with left- and right-moving excitations on [—£, £] is folded into a two-channel c = 4 theory (b) on 
[0,£]. By analytic continuation, this may be written as a chiral c = 4 system on [—£,£] with only 
left-moving excitations. 

FIG. 4. In Euclidean space-time, the finite-size theory (a) on the strip [0,£] with boundary 
condition A at both ends is conformally mapped to a semi-infinite plane (b) with the same boundary 
condition A applied at the boundary. 

FIG. 5. Finite-temperature geometry with periodic boundary condition in imaginary time r. 
The temperature T is identified as 1/(3 and the circumference v/3 of the cylinder satisfies v/3 -C £■ 
For £ — ► oo, this geometry is related to the zero-temperature geometry in Fig. 4b via the conformal 
map w = (v(3/7r)arctan(z). 
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